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Multiple Antenna Communication with Time Reversal Mirror (MIMO-TRM) 

Pre-Processing 

Abstract 

We have pioneered the study of time reversal (TR) pre-processing in radio wireless 
communication. Our results have shown that TR is beneficial for use in a rich scattering 
environment with many propagation paths, thus positioning the technique as a low 
complexity solution for high data rate communication via large delay spread channels. 
Methods to improve the performance of TR have been studied, and comparisons with 
other precoding schemes have been made. In particular, a TR system prototype has been 
designed, and is being built for carrying out fUture experiments. 

Introduction 

Communication via a large delay spread channel at high data rate has proven to be a 
challenging task, primarily due to the substantial amount of intersymbol interference (ISI) 
incurred by the long channel impulse response (CIR). Time reversal (TR) is a simple pre¬ 
processing technique that makes such communication feasible without the need for any 
complex receivers. In a TR communication system, the receiver first sends a training 
sequence to the transmitter for channel estimation. The transmitter convolves the message 
with the time-reversed CIR, and sends it to the receiver. Three main benefits result from 
this simple transmission scheme. 

• Spatial focusing - The spatial profile of the received signal’s power peaks at the 
receiver and decays rapidly away from the receiver. This means that co-channel 
interference in a multi-cell system is low, and the probability of interception by an 
unintended receiver is small. As a result, efficient frequency reuse can be realized, 
and secure communication is attained. 

• Temporal focusing - Analogous to traditional matched filtering at the receiver, 
the power of all multipaths in the CIR can be collected at one particular time 
instant. Thus, the receiver can detect the transmitted message at this time instant 
with maximum signal-to-noise ratio (SNR). 

• Channel hardening - The time variation of the effective CIR seen by the receiver 
is significantly less than that of the original channel. In other words, a high 
frequency diversity gain is achieved, and the adverse effect of fading is reduced. 

Although TR has been studied extensively, and applied successfully, in the ultrasound 
and underwater acoustic arenas, especially in a rich scattering environment with ample 
multipaths, its suitability for wireless radio has yet to be determined. The core of our 
research activities was to investigate the feasibility of applying TR to high data rate 
wireless radio communication. The research results are summarized in the following 
section. For detailed explanations, please refer to the attached publications. 


Summary of Research Results 

1. M. Emami, J. Hansen, A. Kim, G. Papanicolaou, A. Paulraj, D. Cheung, and C. 
Prettie, “Predicted Time Reversal Performance in Wireless Communications Using 
Channel Measurements,” to appear in IEEE Communications Letters. 

This work represents the first effort in estimating the performance of TR in wireless radio 
using real channel measurements. Data were taken in an indoor environment in the 2 to 8 
GHz frequency band by Intel Corporation. Our calculations show that the 3 main benefits 
of TR, namely spatial focusing, temporal focusing and channel hardening, can all be 
obtained. These results support the theory that TR is an equally attractive technology for 
wireless radio communication as for ultrasound and underwater acoustic communication. 

2. J. Hansen, D. Baum, A. Paulraj, “Design Approach for a Time Reversal Test Bed for 
Radio Channels,” Special Session on MIMO Prototyping, 12th European Signal 
Processing Conference, Sept. 2004. 

This work describes the design of a prototype system for carrying out TR experiments in 
wireless radio. The design is heavily based on an off-the-shelf multiple-input-multiple- 
output (MIMO) channel sounder, in order to keep implementation and development costs 
low. A company has tentatively agreed to build this prototype, pending final approval of 
contract by Stanford University and the company. The availability of this prototype will 
allow us to demonstrate that the predicted benefits of TR can actually be attained in 
practical environments. 

3. T. Strohmer, M. Emami, J. Hansen, G. Papanicolaou, A. Paulraj, “Application of 
Time Reversal with MMSE Equalizer to UWB Communications,” Global 
Telecommunications Conference, Dec. 2004. 

This work examines the performance gain of adding a minimum mean square error 
(MMSE) linear equalizer at the receiver of a TR system. Performance improvement is 
obtained because ISI still exists in the received signal, even though the signal is already 
time compressed. Moreover, for fixed equalizer complexity, the concept of rate-backoff 
is introduced. The idea is to reduce the transmission rate, and hence the ISI, to a point 
such that a certain performance level is maintained. 

4. M. Emami, M. Vu, J. Hansen, A. Paulraj, G. Papanicolaou, “Matched Filtering with 
Rate Back-off for Low Complexity Communications in Very Large Delay Spread 
Channels,” submitted to Asilomar Conference on Signals, Systems, and Computers, 
Nov. 2004. (Extended Abstract) 

This work derives the optimum MMSE linear prefilter for an ISI channel, subject to the 
constraint of a simple, one-tap equalizer at the receiver. The performance of this 
transmission strategy is compared to that of TR. 



5. C. Oestges, J. Hansen, M. Emami, A. Paulraj, G. Papanicolaou, “Time Reversal 
Techniques for Broadband Wireless Communications,” European Microwave Week, 
Oct. 2004. (Invited Paper) 

This paper gives an overview of the TR technology as applied to wireless radio 
communication. It contains both a tutorial introduction and recent research results, 
including results discussed in the above contributions. 

Conclusions and Future Work 

We have initiated the investigation on applying TR to wireless radio communication. Our 
results have shown that TR is a promising technology to pursue, with the capability of 
achieving high data rate at low complexity under hostile environment with large amount 
of 1ST Future research will continue to examine TR in details. Questions such as how 
inaccurate channel information affects performance, and how spatial focusing can be 
quantitatively measured, will be addressed. Comparisons with nonlinear precoding 
techniques will also be made. Moreover, we expect to perform TR experiments in 
practical environments once our prototype system is available. On a broader perspective, 
we seek to better understand the fundamental tradeoff between spectral efficiency, power 
consumption, performance and complexity when communicating via large delay spread 
channels at high data rate. 
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ABSTRACT 

In time reversal, a message signal to be transmitted is first 
convolved with the time-reversed channel impulse response 
and then sent to the receiver. The “time-reversed” waves 
which propagate into the channel retrace their former paths; 
they eventually lead to a focus of power in space and time at 
the receiver. Time reversal systems can use antenna arrays 
in order to enhance the fopusing, but they do not necessarily 
need to. The channel’s degrees of freedom are rather excited 
by using broadband signals, and the focusing can be achieved 
with single-antenna links. 

To the authors’ best knowledge, there has been no 
demonstration of time reversal in wireless radio yet. This pa¬ 
per describes the broad design principles of a system that is 
designed to demonstrate time reversal in wireless radio chan¬ 
nels. The system block diagram, some critical specifications 
and sources of error are described. Experimental data will be 
presented, if available in time, at the workshop. 

1. INTRODUCTION 

In a practical system employing time reversal (TR), an in- 
tented receiver sends a training sequence to an intended 
transmitter which is equipped with one or potentially more 
than one transmit antenna. The transmitter time-reverses the 
estimated channel impulse response (CIR) and convolves this 
signal with any message signal that is then sent to the re¬ 
ceiver. 

In quasi-static, reciprocal channels this simple precod¬ 
ing scheme yields a concentration of power at only the in- 
tended receiver at a particular time. The spatial and temporal 
focusing that can be achieved by TR has been successfully 
demonstrated in ultra-sound by Fink [1,2] and in underwater 
acoustics [2,3, 4, 5]. 

In wireless communications, no demonstration is known 
to the authors. In this paper, we propose to convert a com- 
merically available broadband MIMO channel sounder into 
a device that can demonstrate the leverage of TR in low 
Doppler (quasi time-invariant) channels with about one mi¬ 
crosecond delay spread. Our strategy is to keep the upgrade 
as simple as possible. We rely on the MIMO sounder’s capa¬ 
bility for channel estimation, and add the features to transfer 
the CIR estimates back to the transmitter via an external link 
and to then transmit a precoded message. Hence, we replace 
the pair of transceivers usually considered in TR experiments 
by a single transmitter and a single receiver which are con¬ 
nected by some feedback loop. The TR sounder is designed 
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such that this feedback loop can be chosen from a variety of 
publicly available, wired or wireless, links. 

The rationals for this strategy are twofold. For one, the 
modification of this type of existing and verified system min¬ 
imizes the amount of additional hardware effort needed and 
the risk of system performance failure. Secondly, the usage 
of the same link for channel estimation and TR eliminates 
any problems with the reciprocity assumptions that possibly 
occur due to differences in transmit and receiver RF chains 
of a single transceiver. 

The remainder of the paper is organized as follows. In the 
next section, we give a brief introduction to TR. In the third 
section, we describe the system’s design principles. We pro¬ 
ceed to investigate critial specifications and sources of error 
in the fourth, and conclude in the fifth section. 

2. TIME REVERSAL 

The transmitter uses the time-reversed complex conjugate 
of the CIR as the transmit prefilter. We denote the CIR by 
h(rQ, t), where ro is the receiver location and T is the delay 
variable. Applying /C(r 0 ,-T) as the prefilter, the effective 
channel to any location r is thus given by the time-reversed 
field 

s(r,T) =h*(r 0 ,—T)®h(r,T) (1) 

where denotes convolution. 

A convolution with a time-reversed signal is equivalent 
to a correlation. We see from (1) that the focusing relies on 
the decorrelation of CIRs in the delay and the spatial domain. 
A CIR can be considered a random code sequence that is as¬ 
signed to a transmit-receive pair. The auto- and crosscorrela¬ 
tion properties of this sequence are given by nature. In con¬ 
trast to CDMA, they come without any additional bandwidth 
spreading. The decorrelation of die CIRs is obviously best in 
a rich scattering environment. In addition, high bandwidths 
unleash the channel’s degrees of freedom and aid decorrela¬ 
tion. Physically, one can consider the ubiquitous locations of 
the scatterers to form a large virtual aperture which enables 
focusing even with a single antenna down to the diffraction 
limit. A simple but effective measure for the focusing capa¬ 
bility of a TR channel is the delay-spread bandwidth product 
which roughly gives the number of taps that the CIRs has. 
For an initial evaluation of the benefits of TR we refer to [6]. 

3. SOUNDER ARCHITECTURE 

The commercially available MIMO sounding unit consists of 
the following parts. The transmitter contains a signal gener¬ 
ating unit (SGTx) and a transmitter RF module (RFTx), the 
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Figure 1: A schematic of the TR sounder. 


receiver contains a receiver RF module (RFRx) and a data 
storage unit (DSRx). The conversion of this sounder into a 
TR sounder is done by adding components which we call 
TR signal processing unit (TRSPU) and feedback loop (FL). 
Fig. 1 shows the architecture of the TR sounder. 

The proposed TR scheme operates in three steps; the 
main function of the TRSPU is the flexible definition and 
coordination of these steps. In a first step, called channel es¬ 
timation step (CES), the transmitter transmits a training se¬ 
quence that the receiver uses to estimate the channel. This 
estimate is time-reversed in the second step, convolved with 
the transmit sequences, and sent back to die transmitter via 
the FL. Then, in a third step called time reversal step (TRS), 
the transmitter transmits the time-reversed sequence into the 
channel; the receiver now sees a signal that is focused in 
space and time. The time focusing can be demonstrated by 
the compressed C1R at the intended receive antenna; the spa¬ 
tial focusing appears because apart from this one antenna, all 
other ones will detect signals that look like noise. 

In the next subsections we describe the specifications of 
the TR sounder in more detail. 


3.1 MIMO Sounding Unit and Antennas 

3.1.1 Bandwidth , Carrier Frequency, and Output Power 

The maximum bandwidth of the currently available system 
is 120 MHz. For a delay spread in the range of up to 1 fis, 
we can hence expect a time bandwidth product of up to 100. 
Under normal conditions, we expect delay-spreads which are 
a factor of 10 lower than that; hence, the systems works with 
a time bandwidth product of about 10 which still suffices to 
achieve a reasonable benefit with TR. 

The carrier frequency of the system is in the ISM band 
at 2.4 GHz. Unfortunately, the available transmsission band¬ 
width is limited to 83 MHz. According to FCC regulations 
on RF devices [7], the maximum allowed EIRP is 4 W, the 
maximum transmit power is 1 W. The additional 6 dB can be 
achieved by a suitable choice of transmit antenna. 

A second suitable choice would be the ISM band at 5.8 
GHz with an available bandwidth of 150 MHz. However, as 
the free space attenuation is higher and thus reduces trans¬ 
mission range, we decided to perform the demonstration at 
the lower carrier frequency. 

The output power of the system is limited to 1 W. 


3.1.2 Signaling 

The sounding signals, and in particular their power and 
power spectral density distribution, are designed in fre¬ 
quency domain. The signals used for channel estimation 
have constant maximum amplitude over frequency and their 
phase is designed to minimize the crest factor and hence 
maximize average output power. Since all signal processing 
is performed in frequency domain, instead of estimating and 
time-reversing the CIR, the transfer function (TF) is complex 
conjugated in frequency domain. 

3.1.3 Antennas and Switching 

In our system, the antenna elements in an array of the sounder 
are switched. In this case the number of channels is essen¬ 
tially limited by the number of available switches and an¬ 
tennas and the maximum measurement time for which the 
quasi-static channel assumption still holds. For the TR ex¬ 
periments, arrays with 4 and 8 antennas will be employed. 

Since TR requires a large angular spread, omnidirectional 
antennas are most advantageous. However, the range that is 
reqiured to be large in order to obtain a large delay spread 
would benefit from high gain antennas. The compromise be¬ 
tween these two constraints are basestation type transmit an¬ 
tennas (very narrow vertical and wide horizontal beamwidth, 
gain about 12dBi) and slightly directional antennas (very 
wide horizontal and wide vertical beamwidth, gain of 6 to 
10 dB). 

3.1.4 Synchronization 

The synchronization is done by using very stable reference 
clocks at 10 MHz (e.g. Rubidium standards). The clock sig¬ 
nal guarantees longterm stability. In each RF unit, the carrier 
signal is generated independently by low-noise PLLs. 

3.2 Feedback Loop 

In between the CES and the TRS, the estimated channel data 
has to be transferred from the receiver to the transmitter. The 
transferral time must be much shorter than the coherence 
time of the channel. The data rate of the link becomes hence 
a critical parameter. A single complex CIR with about 1000 
taps and 8 bit resolution has a size of about 16 kbit. Assum¬ 
ing a time invariant channel with Doppler frequency of 1 Hz, 
a realistic transmission time should not be longer than 0.2 s; 
consequently, the FL needs a data rate of at least 80 kbit/s per 
CIR that is transferred. 

The following options have been evaluated. 

3.2.1 Microwave Link 

A microwave link is advantageous in terms of data rate, 
which is in the order of Mbit/s, and in terms of flexibility. Mi¬ 
crowave links can be established whereever the experiment is 
to be performed. However, TR experiments are likely to be 
conducted in environments with large delay spreads. Here, 
the TR sounder will operate at maximum EIRP allowed by 
the FCC regulations. This limit holds for the microwave link, 
as well. In areas where the SNR of the TR experiment is low, 
that of the microwave FL is low, too. The link quality can 
be improved by using directional antennas mounted on high 
poles, but this requires difficult adjustments and large efforts 
in the preparation of the experiment 








3.2.2 Cellular Systems 

Currently available cellular systems such as GPRS offer great 
flexibility without the need for heavy experimental equip¬ 
ment. Their main drawback is their comparably low data rate 
(order of 100 kbit/s for GPRS) and, in particular, the unreli¬ 
able transmission delay. Since the data rate offered does not 
provide much overhead compared to the data rate required, 
the application of current cellular services for the FL appears 
risky, in particular if the Doppler spread of the channel may 
at times exceed 1 Hz. 

3.2.3 Wired Network 

The wired network offers tremendous data rates; the main 
drawback is the inflexibility of its access. However, most TR 
experiments will be conducted in urban areas, where wired 
access to the internet can be expected to be in range. A wire¬ 
less network system (WLAN) such as one based on the IEEE 
802. lib standard still offers data rates in the order of Mbits 
and can be employed to bridge the gap between the experi¬ 
mental device to the next wired network connection. 

The final decision about the best link depends on the 
availability of particular services, on local network traffic, 
and the Doppler spread of the channel. The sounder itself is 
equipped with RJ45 connectors and can be hooked up to any 
of the above-mentioned systems. 

3.3 The TRSPU 

The role of the TRSPU is to coordinate the three steps on 
which the operation of the sounder is based. Since it also 
performs the nessecary signal processing, it is located in be¬ 
tween the RF and the data storage unit of the receiver. 

3.3.1 The Different TR Modes 

The tasks of a TR sounder are at least threefold: 1) It must 
demonstrate spatial and temporal focusing when TR exper¬ 
iments are performed with a) a single, or b) multiple trans¬ 
mit antennas. 2) It must be flexible to allow further research 
about benefits of TR, i.e. for multiuser (MU) communica¬ 
tions. 3) It should be backwards compatible. Since it is based 
on commerical MIMO technology, MIMO channel measure¬ 
ments should still be possible with a TR sounder. It turns out 
that 4 different TR modes evolve from these requirements: 

1. MIMO Measurements: 

For MIMO measurements, the TRS is omitted. The 
sounder performs a predefined number of channel esti¬ 
mations and no TR transmissions. 

2. SISO-TR: 

In the CES, the channel estimation is performed times 
with a single transmit and a single receive antenna. The 
received TFs can be averaged to obtain a better estima¬ 
tion of the channel. The TRSPU complex conjugates this 
TF, multiplies it with the Fourier transform of the used 
transmit sequences, and sends it via the FL to the trans¬ 
mitter. The transmitter loads this sequence and transmits 
its Fourier-inverse into the channel. Now, in the TRS, the 
TRSPU provides a second switching table and the RFRx 
switches through the antenna array. At the one receive 
antenna that is used in the CES, the detected CIR is com¬ 
pressed. At all other antennas, a noise-like signal ap¬ 
pears. After «2 transmissions, the CES is repeated. This 
mode is the simplest and best suited to demonstrate TR. 


SGTx 


TRSPU 



Figure 2: The lower layer of the TRSPU, synchronized with 
the SGTx. 


3. MU-SIMO-TR: 

In this mode, MU communication with TR can be 
demonstrated. The CES is performed with a single trans¬ 
mit antenna and an antenna array at the receiver. Concep¬ 
tually, it works similar to the SISO-TR mode, except that 
the TFs of all receive antennas are estimated. The TF es¬ 
timated at each antenna is multiplied with a message that 
is to be sent to this antenna; the sum of all these products 
is fed into the FL and retransmitted in the TR-mode. The 
channel acts now as a filter; each receive antenna receives 
only its intended message, and all the other messages as 
noise. 

4. MISO-TR: 

In the original underwater acoustic applications, TR is 
performed with a transmit array. Each transmit antenna 
focuses its signal on the same receive antenna, so that 
the focusing is enhanced proportional to the number of 
transmit antennas. Since in our application, the arrays 
are switched, the superposition of the signals can only 
be done by software. Also, the switching requires short¬ 
term phase stability at the transmitter. 

This 4th mode is not only the most demanding with re¬ 
spect to synchronization, but also with respect to the data 
rate on the FL, since as many CIRs must be transferred 
as there are transmit antennas. 

3.3.2 Double Layered Structure of the TRSPU 

The TRSPU is implemented in a double-layered structure. 
The lower layer runs synchronized with the SGTx and is as 
such displayed in Fig. 2. The SGTx receives a sequence from 
the FL, transforms it into time domain and hands it over to 
the RFTx. The TRSPU receives the demodulated signal from 
the RFRx and computes the conjugate of the received CER 
via an FFT and division by the sequence used for channel 
estimation. This layer is run in all modes of the sounder, be 
it MIMO channel estimation, or one of the two stages of the 
different TR modes. 

The upper layer bascially determines the TR mode, as 
described in Subsection 3.3.1. Each mode is defined by par¬ 
ticular numbers n\ and «2 for the duration of the CES and 
the TRS, and the one (for the MIMO mode) or the two (for 
all other TR modes) corresponding switching tables that are 
provided to the RFRx. The upper layer is software-defined. 
This allows any user to make changes in the setup; in partic¬ 
ular, the TR modes can be embedded in more sophisticated 
high-level transmission schemes, e.g., schemes that adapt the 
numbers n\ and «2 to particular channel conditions. 

An example of the SISO-TR mode is shown in Fig. 3. For 



















1) The channel estimation loop is run n 1 times 



2) At the beginning of the transmission 
mode, the CSi is sent to the Tx _ 

I .J 

3) In the transmission mode, predefined sequences ire 
transmitted 


Figure 3: The SISO mode as triggered by the upper layer of 
the TRSPU. 


a given transmit antenna Tx(0) and an intended receive an¬ 
tenna Rx(0), the TF H((0, (0,0) ) is estimated in the first step. 
In the second step, it is multiplied with die Fourier transform 
of a transmit sequence and transmitted to the FL. In the third 
step, the sounder transmits these combined sequences until it 
again enters the first step. 

4. CRITICAL SPECIFICATIONS AND SOURCES 

OF ERROR 


43 Phase Noise 

TR requires exact phase information at the transmitter. Inac¬ 
curacies in the phases occur since a) the channel is not per¬ 
fectly time variant, b) the LO signals derived by PLLs from 
the 10 MHz reference have phase noise. 

The time variance can be overvome if the data rate of the 
FL is large enough (Subsection 3.2 and 4.2). 

Phase noise destroys the coherent superposition of sig¬ 
nals transmitted from multiple antennas. Since it enters both 
at the transmitter and the receiver, it will reduce the received 
power by a factor of cos(a<£i 4-a <fc), where a<& - / = 1,2 are 
random variables that describe the phase fluctuations of the 
transmitter’s and the receiver’s oscillators. For phase errors 
less than 15°, this yields a performance loss of about 10 % in 
SNR, which is tolerable. For switched antenna sounders this 
is less of a problem with respect to TR. 

5* CONCLUSIONS 

This paper introduces time reversal (TR) as a promising 
transmission technique for wireless communications. We 
propose the design of a prototype device that will allow the 
demonstration of the benefits of TR. This prototype is heavily 
based on commercially available systems and components, 
implementation and development costs are kept at a mini¬ 
mum. The core of the device is a commercially available 
MIMO broadband channel sounder. It can be connected to 
various types of publicly available networks which serve as 
a feedback loop. Its key component is the TR signal pro¬ 
cessing unit, which can operate in several modes and offers 
flexible handling of the entire device. 
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Abstract 

We propose to apply a technique called time-reversal to UWB communications. In time-reversal a signal is precoded such 
that it focuses both in time and in space at a particular receiver. Spatial focusing reduces interference to other co-existing systems. 
Due to temporal focusing, the received power is concentrated within a few taps and the task of equalizer design becomes much 
simpler than without focusing. Furthermore, temporal focusing allows a large increase in transmission rate compared to schemes 
that let the impulse response ring out before the next symbol is sent. Our paper introduces time-reversal, investigates the benefit 
of temporal focusing, and examines the performance of an MMSE-TR equalizer in an UWB channel. 

Index Terms 

Time Reversal, MMSE Equalization, Ultra-Wideband, Intersymbol Interference 


Hansen’s work was supported by the Deutsche Forschungsgemeinschaft. This work is partially supported by grants AFOSR F49620-01-1-0465, ONR 
N00014-02-1 -0088, NSF DMS grant 0208568 and by NSF Contracts CCR-0241921, 0241919, ONR Contracts N00014-03-1-0766, N00014-03-1-0906 and 
Intel Corp. 





1 




i 


I. Introduction 

Ultra-wideband (UWB) has become a suitable candidate 
for high-data rate, short range communications. Due to the 
large operation bandwidth, the resolution in delay domain is 
extraordinary, so that even in a dense-scattering environment, 
the probability of fading is low. Since for the generation 
of these ultra-short pulses, no filters are required, it was 
considered that UWB devices are cheap to produce. 

Recently, however, several drawbacks have been noted. 
Even though UWB transmitters irradiate very low power per 
bandwidth, they can potentially interfere with many systems 
that have frequencies assigned in the low GHz range [1]. 
Furthermore, multipath fading may be low, but in order to 
catch about half of the energy that is distributed in the 
entire impulse response, RAKE receivers with at least 20, but 
potentially many more taps must be constructed [2], [3]; for 
handsets, such a design is not low-cost. 

We claim that both these drawbacks can be overcome 
if UWB is combined with a transmission scheme which is 
called time reversal (TR). This scheme has its origin in wide¬ 
band transmission in under-water acoustics [4], [5], [6], [7] 
and ultra-sound [8], [7] and has recently attracted attention 
of wireless communications engineers [9]. In TR, the time- 
reversed channel impulse response (CIR) of any transmit- 
receive link is taken as a prefilter at the transmitter. If such 
a time-reversed sequence is irradiated into the channel, its 
components retrace their former paths and lead to a focus 
of power at the intended receiver at some particular time 
instant. For UWB communications, this strategy has several 
advantages. Focusing power means that the receiver needs only 
very few taps to capture a significant amount of the power 
in the channel. Essentially, the complex taks of estimating a 
large number of channel taps is removed from the receiver, 
which may be a low-cost handset, to the transmitter. In 
addition, the rate of the system can be increased. Whereas with 
proposed pulse-position modulation schemes the time between 
the repetition of two pulses is chosen to be sufficiently larger 
than the length of the CIR [10], focusing of power means 
that the repetition rate can be increased. At last, the spatial 
focusing that comes along with TR yields lower interference 
with other communication systems. 

In this paper, we apply TR to a single user downlink 
scenario of an UWB channel. In order to assess the usefulness 
of TR, we define a time-compression factor which measures 
the fraction of the energy in the non-focused components 
of the channel versus that in the main, focused peak. For a 
channel with iid complex zero mean circularly symmetric taps, 
50% percent of the entire CIR’s energy can be captured with a 
single tap receiver. In real channels, this ratio is slightly lower, 
which we show with the aid of experimental data from UWB 
measurements. 

Since the scattered components are usually widely dis¬ 
tributed over the delay axis, a receiver that wants to capture 
more than 50% of the signal energy in the effective CIR must 
use some equalization technique. We propose an TR-MMSE 
equalizer for this purpose. We investigate the performance of 
this equalizer on real channel data and show that it can operate 


at sampling rates which are much higher than the inverse of 
the length of the CIR of the channel; no ringing out of the 
impulse response is required. 

This paper is organized as follows: in the second section, we 
introduce TR and define and investigate the time compression 
factor. Also, the signal model is described and the principles 
of the TR-MMSE equalizer are outlined. In the third section, 
we first investigate the time compression factor with the aid 
of real data. Then, we examine the performance of the TR- 
MMSE equalizer in term of bit-error-rates (BER). In particular, 
we vary the number of taps and the sampling rate. We conclude 
in Section IV. 

II. Theory 

A. Time Reversal 

The transmitter uses the time reversed complex conjugate 
of the CIR as the transmit prefilter. We denote this CIR by 
/t(ro, t), where ro is the receiver location and r is the delay 
variable. Applying the complex conjugate /i*(ro,— r) as the 
prefilter, the effective channel to any location r is thus given 
by the time reversed field 

g(r,T) =-^ h *{r 0 ,-T)®h(T,T) (1) 

where g> denotes convolution and the factor normalizes 
the transmit power with the square root of tire channel’s 
energy, E H = / |/i(r 0 ,r)| 2 dr. 

A convolution with a time-reversed signal is equivalent to 
a correlation. We see from (1) that the focusing relies on the 
decorrelation of CIRs in the delay and the spatial domain. In 
a rich scattering environment, a CIR can be considered as a 
random code sequence that is assigned to any transmit-receive 
pair. The auto- and crosscorrelation properties of this sequence 
are given by nature; unlike in CDMA, they come without 
any additional bandwidth spreading. In UWB, these sequences 
are particularly long. The benefits of time compression are 
crucial, and the shortening of the CIR reduces the complexity 
of equalization at the receiver. 

The decorrelation of the CIRs is obviously best in a 
rich scattering environment. Physically, one can consider the 
ubiquitous locations of the scatterers to form a large virtual 
aperture which enables focusing even with a single antenna 
down to die diffraction limit. A simple but effective measure 
for the focusing capability of a TR channel is the delay spread- 
bandwidth product which roughly gives the number of taps 
that the CIRs has. For an initial evaluation of the benefits of 
TR we refer to [9]. 

B. Temporal Focusing 

In a multipath environment, the channel’s energy Eh is 
spread over delay. It can be decomposed into a term Eq 
which describes the fraction of the energy in the direct path, 
and a term E$, describing that of the multipaths. In a strong 
multipath environment such as given for UWB, we have the 
ratio Es/Eq 1. If the channel is slowly time-variant, the 
quantities E$, Eq, and Eh are random variables. Their mean 
is denoted by (.). 
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In TR, the CIR is compressed and a temporal focus of 
energy is visible in the center of the compressed CIR. In order 
to characterize the amount of the temporal focusing, we define 
the temporal compression ratio as 

7Tfl = (E T s R )l(1E[™) = - 1 (2) 

where E$ R is the energy in the main peak of the received 
impulse response, Eg R the one in the tails, and E R R is the 
sum of the two. 

We can compute 7 tr in terms of either the CIR h(r) or 
the transfer function From (1) we see that the transfer 

function H tr (cj) of the time-reversed channel is the squared 
absolute value of that of The discrete CIR h TR [l], with 
a total of L non-zero taps, is 

h TR [l] = + J |//(w)| 2 exp(jwl)du>. (3) 

From this we compute the energy of the Oth tap of the time- 
reversed channel, 

(e™) = (i^toii 2 ) - a£ f 

= {( E s + E 0 f) = (E 2 H ) (4) 


and the energy of the entire channel, 

(Ejj R ) = ($2 Ml ® V[~l}\ 2 ) = I# HI 4 <M (5) 


where the frequency domain representation on the right hand 
side of the equation follows from Parseval’s theorem. Whether 
a time domain or a frequency domain representation is more 
desirable depends on in which domain information about the 
channel is available. 

If we evaluate expression (5) further in time domain for L 
identical and independently distributed taps with equal power, 
we can, derive the time compression factor as 

04 + 2 (L - l)/% 


7 th = 


04 + (L - 1)01 


- 1 . 


( 6 ) 


Here, 0 2 and 0\ denote the second and the fourth moment of 
the amplitude in each tap, respectively. For L — 1, one has 
7 TR = 0; with increasing L it converges to 1. Hence, in the 
time-reversal channel, the ratio of the energy of the zeroth tap 
to the sum of the energy of all the other taps converges to a 
fixed value. 

If h[l] is a complex Gaussian random variable with mean 
0 and variance {Eq + Es) = {Eh), we have 0 2 = {Eh) and 
04 — 2 {Eh) 2 - Hence, 


ITR 


2L{E h ) 2 _2 

2{E h ) 2 + {L-1){E h ) 2 1 + i 



For large L, the power in the tails of the impulse response 
becomes about equal the power in the peak. Hence, even 
though ISI can be significantly suppressed, any system without 
equalization operates at best at 0 dB signal to ISI plus noise 
ratio. For small L, the compression works in the mean better 
than for a large one, even though the difference quickly levels 



Fig. 1. Schematic of TR impulse h tr and truncated TR CIR hj, R . In this 
toy example L — 29, L s — 9 and the number of effective taps (marked as 
fi Ued circles) is 5. 


out. But the probability that a channel realization occurs that 
does not focus at all (i.e., all channel taps are real and have 
the same magnitude) is much larger than for long sequences. 

C. Signal Model and TR-MMSE Equalizer 

The transmitter modulates symbols Sk using a pulse shaping 
filter 0(r) of bandwidth B — l/T s . It then performs TR on 
the pulse stream. Thus: 

x(t) = ^ s k (j){t — kMT s ) <g) h*(—t) (8) 

k 

where h{r) is the (infinite bandwidth) CIR. After matched 
filtering the received signal can be written as: 

y{t) = — kMT 3 ) ®h*(—t)®h(t) <&4>*(-t) + n(t) 

k 

(9) 

The receiver then samples y(t) at IMT S : 

y[l] = y(lMT s ) = ]T ~k}+ n[l], (10) 

k 

where limit) = h(t) <g> h*(—t), h{t) = h(t) <g> <p{t), and 
hm[k] = hm{kMT s ). The sampled received signal is thus 
the convolution of the symbols with a downsampled-by-M 
version of 

Due to the temporal focusing described above the equalizer 
at the receiver does not have to use the quite long (time- 
reversed) CIR. Instead the equalizer uses a “shortened” TR 
CIR by keeping only a few taps which capture most of the 
energy of hm and still achieve very good performance. More 
precisely, we construct a shortened TR CIR h^ R of length L s 
as illustrated in Fig. 1 by defining 

*™i<]» {j tr[!i on 

for |/| < L s /2 and L s — max{Z : \Iitr[1]\ > 6 }. Here 6 is 
an a priori chosen tolerance level that depends on the target 
SNR. For convenience we will rename the indexing of h RR 
such that hjiji = {^tr [ 0 ]> • • •, h^ R [L s — 1 ]]*. 
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The symbols are transmitted one block at a time with a 
guard period of length L s . The output of the channel is: 

y[m] = (x 0 /ir;?)[ra] + n R {m\ 

Here x[m\ is the input signal m ~ 0,..., N — 1 and y[m] is 
the output signal m = 0,..., N -f L — 2. In matrix notation 
this can be written as Toeplitz-type system 


y = HrflX + n R 


where y — [t/[0]... y[N -f L - 2]] T and x = [z[0]... x[AT - 
1]] T . Using hj> R instead of Iitr for the MMSE equalizer we 
replace H tr by the matrix Hf R of size (N + L s — 1) x N 
given by 

f h a TR [0] 0 ... 0 1 


m 

* 

0 h^ R [L 8 — 1 ] _ 

and the matrix equation for the MMSE estimator becomes 

X = (W*)*H^ + (H^)*y. 

Assuming the transmitted symbols are from a BPSK alphabet, 
we can now decode x by looking at the sign of each element. 

We note that is a sparse Toeplitz matrix, thus x can 
be computed efficiently by combining sparse matrix techniques 
with fast Toeplitz solvers [11]. 


— 


^TrI^s 

0 


0 


III. Simulations 

A, Measurement Data and Temporal Focusing 

Measurements were conducted with a network analyzer by 
Intel Corp. at off-peak hours to ensure channel stationarity. 
The environment is an office space (40m x 60m) with many 
cubicles; measurements were conducted at several locations. 
They span a bandwidth of 2-8 GHz with 3.75 MHz frequency 
resolution. Antennas are vertically polarized. The measure¬ 
ments are corrected to compensate for the system components 
(including cable, gain stages, and antennas). The height of 
the transmit antenna is about 2.5m and that of the receive 
antenna is 1m above the floor. As a sample result, Fig. 2 
shows the CIR (above) and the time-reversed CIR (below). 
The compression is clearly visible. Still, this figure is not 
suitable to tell the ratio between the power in the peak of 
the time-reversed impulse response and that in the tails. In 
order to demonstrate this ratio, we plotted in Fig. 3 the time 
compression ratio j tr for a CIR length between 1 and 500 
taps, where the entire length of the CIR is fixed at l/3.75MHz. 
The curve was actually computed in frequency domain; in 
order to suppress fluctuations due to the randomness of the 
samples, a sliding window was implemented, so that ^tr was 
actually averaged over all available frequency samples. 



Fig. 2. Impulse Response of the channel (above) and of the time-reversed 
channel (below) 



Fig. 3. In the TR channel, the ratio between the eneregy in die main tap and 
the energy in the tail of the CIR converges for independent taps to a fi xed 
ratio. 


For a small number of taps, the curve coincides very well 
with the theoretical one. At about 200 taps, the maximum 
theoretical value of 1 is exceeded; beyond that, correlation 
between the taps starts to play a role and the power in the 
tails of the CIR will be stronger than that in the main peak. 


B. Performance of the Equalizer 

We used the single channel realization whose magnitude 
response at full bandwidth, B, is shown in Fig. 2. Three 
cases were considered. Case I has M = 8 (see Eqn. 10) and 
bandwidth B. Case II has M — 4 and bandwidth B/2. Finally 
case HI has M — 2 and bandwidth B/4. The constellation 
was BPSK and no channel coding was used. Thus all the 
scenarios have the same data rate and all channels have 250 
taps. Fig. 4 shows the bit error rate (BER) curves for these 
scenarios when the receiver only estimates the 20 strongest 
taps. When M = 2 the curve floors very rapidly since not a 
very significant portion of the energy of the CIR is captured by 
the largest 20 taps. When M = 4, the ISI energy is due only to 
every fourth tap of the TR channel. Thus more of the energy 
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Fig. 4. BER vs SNR when the receiver only estimates the 20 largest taps 
of the channel. Bandwidth expansion factor of 8 (solid line), 4 (dotted line), 
and 2 (dashed line) 



Fig. 5. BER vs SNR when the receiver only estimates the 10 largest taps 
of the channel. Bandwidth expansion factor of 8 (solid line), 4 (dotted line), 
and 2 (dashed line) 


Fig. 6. BER vs SNR when the receiver only estimates 20 largest taps of the 
channel - another channel realization. Bandwidth expansion factor of 8 (solid 
line), 4 (dotted line) and 2 (dashed line). 



Fig. 7. BER vs number of strongest taps estimated at the receiver. SNR = 
15 dB 


of the effective CIR is captured by the strongest 20 taps. This 
effect improves significantly in the case where M = 8. Note 
that the case M — 8 uses more bandwidth. However, the BER 
floors at a reasonable SNR of 15 dB for UWB. Fig. 5 shows 
the BER curves for the same three cases but when the receiver 
estimates only 10 strongest taps of the CIR. In this scenario 
the error flooring occurs earlier for ail the cases. However, the 
complexity of the receiver is significantly reduced. 

Fig. 6 plots the BER for a different CIR. Here the receiver 
estimates the largest 20 taps of the CIR. Comparison between 
this figure and Fig. 4 shows that the BER performance is very 
similar for these two channel realizations. This is due to the 
large diversity available in this ultra-wideband channel. 

Fig. 7 shows the BER vs number of strongest CIR taps that 
the receiver estimates for equalization at 15 dB SNR. This plot 


shows that the BER reduces approximately exponentially with 
increasing number of taps. Fig. 8 shows the BER vs number 
of strongest CIR taps that the receiver estimates at 15 dB 
SNR when the transmitter does not use TR Figs. 7 and 8 
demonstrate the temporal focusing ability of TR The receiver 
requires to estimate much less taps with TR than without it in 
order to achieve the same BER performance. 

Note that the plots are for the strongest taps. These taps 
may not be next to one another. There may be zeros between 
any two consecutive strongest taps. Thus the total length in 
samples of the estimated channel is in general larger than the 
number of strongest taps estimated. However, the complexity 
of equalization is approximately the same as if the length of 
the estimated channel is equal to the strongest number of taps 
estimated. 
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Fig. 8. BER vs number of strongest taps estimated at the receiver with 
regular transmission. SNR - 15 dB 


IV. Conclusion 

Time reversal significantly reduces the number of taps that 
the receiver needs to estimate in order to achieve a certain 
target BER. If the CER has iid complex Gaussian taps the 
ratio of signal to ISI power approaches unity as the length of 
the CIR increases without bound. We show that this ratio is 
very close to unity when the CIR is more than about 10 taps 
long. 

In order to mitigate the ISI one can fix the bandwidth and 
increase symbol spacing by an integer factor of M . The larger 
M y the less taps of the CIR the receiver needs to estimate in 
order to attain a reasonable BER. However, larger M means 
more wasted bandwidth. With no channel coding and by only 
estimating the 20 strongest taps out of approximately 250 
channel taps, the receiver can reach 10" 3 BER at around 12 dB 
SNR with MMSE equalization and M = 8. In order to achieve 
such a performance without TR, the receiver needs to estimate 
almost all the taps of the CIR. 
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Extended Abstract 

We study the possibility to transmit data over channels with 
large delay spreads under the constraint of a very simple re¬ 
ceiver which has only one tap. Such a scheme is of interest 
when a cost-efficient way to transmit potentially large data rates 
is sought. We investigate the performance of die optimal pre¬ 
filter for this scheme, and compare it to a simplified, so-called 
time-reversal, prefilter which has very low complexity. Time 
reversal filters have in connection with wireless transmission 
over large bandwidths their origin in underwater acoustics and 
ultrasound [1], Their application to wireless communications 
was recently demonstrated on UWB data [2]. At low SNR, the 
optimal prefilter and the TR prefilter are equivalent. At high 
SNR, the time reversal prefilter achieves a performance that is 
independent from the delay spread of the channel and hence the 
same for any bandwidth. 

In this contribution we derive analytically the performance 
of the discussed system. We discuss the impact of multiple 
transmit antennas and of a transmission scheme that uses rate 
back-off, i.e., it transmits at a rate lower than the inverse of die 
bandwidth. If bandwidth is abundant, such as it is the case in 
Ultra-wideband communications [3], our proposed scheme can 
achieve a reliable quality of service at high data rates with a 
single tap receiver. We mention that our setup corresponds to 
the downlink of an access point to a handheld. In principle, 
our analysis holds equivalently for a simple transmitter with a 
receive equalizer with a large number of taps. 


A. Prefiltering for Long Delay Spread Channels 

A block diagram of the system under study is given in Fig. 1. 
The transmit symbols x are fed into equalizers gi and then trans¬ 
mitted through the channels described by their channel impulse 
responses (CIRs) h i} i= 1...., Mt, and Mr is the number of 
tr ansmi t antennas. The noise is denoted by n, and the receiver 
has a single tap whose gain is denoted by / 3 . Since with very 
large delay spread the transmit prefilter is not likely to eliminate 
inter signal interference (ISI) completely, we allow that the sys¬ 
tem does not transmit at full rate (inverse bandwidth), but at a 
integer fraction of this rate, which we call D. The input-output 



Fig. 1. A block diagram of a MISO system with a preequalizer at the trans¬ 
mitter and a s ingle tap receiver. 


relation of this system can be written as 

Mt 

y lD] [k\ = p y'(hi * 9i )[Dl]x[k -l}+ I3n[k] (1) 

mmm mot 4*«mp 

*-l l 

where [k] is the output symbol at time k. The rate back-off 
factor D = 1,.... L is equivalent to upsampling the signal be¬ 
fore transmission by a factor of D, and L is the channel length. 

We derive the optimal prefilter solving the minimization 
problem 

Sp,N = argmaXg.^lil x*/3Hg + Pn[k] - x[k - A] || 2 ] 
s t. E || g*g ||= 1. (2) 

The product Hg is the convolution of the prefilter and the chan¬ 
nel, the latter stacked into a toeplitz matrix, and N is the length 
of the filter. The vector g is given by 

g = [g[0] T ... e\N - i] T f (3) 

and 

g[i] = l«iW 9m t [1}} T i = 0,... ,2V — 1. (4) 

A is the delay of the equalizer and is approximately — 2 —. 
The difference between this formulation and the one where the 
equalizer is placed at the receiver is the transmit power con¬ 
straint. The above problem is convex in g and 8 and can be 
solved analytically using the method of Lagrange multipliers: 

g P,N = i (h h H + H H e A (5) 

P 2 = ejH (h h H + ylj H"e A (6) 


f 
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We furthermore define the matched-filter bound SNR Pmfb as 
PMFB = -~2-• W 

Combining ( 8 ) and (9) we have that the effective SNR of a TR 
communication system with Mt transmit antennas and a rate- 
back off of D is 

MUD = PMFB— -(10) 

1 + Pmfbbip“j 

In our paper we will present a formula for the ratio of the use¬ 
ful power against die ISI power, E[Ptr]/ E[Q^] for a chan¬ 
nel with exponentially decaying power delay profile and delay 
spread cr T . Using this result, we derive the following limiting 
cases: 


Fig. 2. The effective channel as generated by TR (dashed line) and the MMSE 
filter (coot line) 


Eq. (5) defines a filter with an MMSE-type solution with a 
finite number of taps. We notice that the optimized value for (3 
bascially ensures that the power constraint of die filter is valid. 
In the transmit-receiver chain of Fig. 1, (3 has no effect on the 
SNR at the receiver; it ensures, however, that the transmitter 
can take into account the SNR to optimize the preequalizer. 

We introduce a simplified filter which we call TR-prefilter 
and which basically is a matched prefilter; it follows from the 
optimal filter (5) for low SNR, 


g TR 


lim g = h* 

SNR-+0 
s.t. || h ||= 1 , 



or 

9TR,i[k\ ~ i/WH, 

with 7 = YliZi l^iMi 2 ) i ie., Ihe filter equals the 

time-reversed CIR of the channel, and the receiver receives the 
autocorrelation function of the channel. The power in its central 
peak is maximized, but the ISI suppression is worse. 

The CIR of a pre-equalized channel with originally 50 taps 
and constant power delay profile is shown in Fig. 2; for the 
computation of die MMSE filter, p = Ptx/& 2 = 5 dB was 
assumed. The dashed line is obtained from the MMSE filter, the 
continuous one from the TR prefilter. We see that die prefilters 
compress the CIRs significantly around a strong peak. The ISI 
of the impulse response obtained with TR is higher than that of 
the MMSE filter, but its peak is slightly higher, too. For a large 
number of taps, the peak is essentially non-fading which yields 
reliable system performance. 


B. Performance Evaluation of the TR filter 

For performance evaluation, we denote by Po and the 
instant useful power and the ISI power, respectively, and define 
the effective SNR p^ff as 


Peff,D = 


Pt*E[P 0 ] 
Pt x E[QW] + 


E[P 0 ] 

E\QV> l] + i‘ 



Peff,D 

Peff,D 


Pmfb for 0 T —> 0 

DMtpmfb r 
-——— for oy —> oo 

Pmfb + DMt 


(H) 

( 12 ) 


Eq. (11) tells that a TR system is lossless should the delay 
spread of the channel be very small. Eq. (12) can be split into a 
low-SNR and a high-SNR regime as 


peff,D = pmfb for pmfb low. <j t large (13) 
peff,D = DMt for pmfb-, cr T laige. (14) 

Eq. (13) shows that TR is optimal at low SNR; this behavior is 
expected for a transmit matched filter. Numerical results pre¬ 
sented in die final paper will show, however, that the low-SNR 
region ends at about 0 dB which not sufficient for reliable com¬ 
munications. Eq. (14) demonstrates that in an environment with 
a large delay spread a T , the performance of a TR system at 
high input SNR saturates at a value independent from the delay 
spread of the channel. The performance increases with 3 dB per 
added transmit antenna, and with another 3 dB for each rate- 
back off of a factor 2 from the maximum available bandwidth. 
We will also include simulations that show that the central tap 
is non-fading. 

If bandwidth is abundant, we can hence design a commu¬ 
nication system with a very simple receiver that has only one 
tap and does not need to perform channel estimation. This tap 
is non-fading and the effective SNR is not a function of the 
bandwidth of the system. Varation of transmission rate or us¬ 
ing some particular number of transmit antennas guarantees any 
reasonable quality of service. In our final paper we will discuss 
the derivations and results presented here more thoroughly, give 
some numerical examples which highlight some of the analyti¬ 
cal results, and point out future applications of the investigated 
system. 
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Abstract —Channels with large delay spreads, which typi¬ 
cally arise when large band widths are used, are considered as 
disadvantageous for wireless communications. However, despite 
their challenges, these channels also offer opportunities. The 
application of time reversal techniques may be one possibility 
to exploit these advantages and to make communications in 
channels with large delay spreads more feasible. This paper 
describes advantages and draw backs of time reversal techniques 
in broadband wireless communications, and outlines concepts, 
goals and potential improvements. The properties of time re¬ 
versal In various simulated and measured wireless channels are 
demonstrated, and system design is discussed. 

I. Introduction 

While there is little doubt that future communication sys¬ 
tems will have to accomodate larger data rates, the natu¬ 
ral solution consisting of increasing the bandwidth brings 
a number of open challenges. In this matter, time reversal 
(TR) techniques could play a role in the design of broadband 
wireless communications. In TR, channel state information 
(CSI) at one or possibly more transmit antennas is used to 
precode transmitted symbols with the time reversed version 
of their respective channel impulse responses (CIRs). These 
“time reversed” waves propagate in the channel, retrace their 
former paths and eventually lead to a focus of power in space 
and time at the receiver [l]-{3]. In rich scattering media, this 
space-time compression can be very strong, and at the focal 
point, the effective channel obtained through TR is hardened. 

Spatial focusing [4] means that the spatial profile of the 
power peaks at the intended receiver and decays rapidly away 
from the receiver. Temporal focusing means that the CIR at the 
receiver has a very short effective length. Channel hardening 
is the phenomenon of tightening of the distribution function 
of the effective CIR’s power. 

Using TR for wireless communications in channels with 
large delay spreads presents therefore three merits. 

« The temporal focusing compresses the channel’s energy 
into one single peak with very low sidelobes. The task of 
equalizing the channel is hence significantly simplified. 

♦ The spatial focusing can result in much lower intercept 
probability and co-channel interference in a multi-cell 
system. By selectively focusing the energy in both space 
and time at a target point, TR ensures that intercept 


receivers will have difficulty detecting or decoding the 
intended signal. 

• Through the channel hardening, TR provides high diver¬ 
sity gain. The statistics of the time reversed channel are 
different from the actual channel. Specifically the time 
reversed channel has a much smaller variance than the 
physical channel itself. 

Time reversal is well known in acoustics and has lead to re¬ 
markable applications in underwater sound [5]—[13] and ultra¬ 
sound [14]—[16]. The extension of TR techniques to wireless 
communications has yet to be precisely investigated, although 
the idea of exploiting scattering has recently been proposed 
in [17]—[20]. From a communication engineer’s viewpoint, TR 
shifts the system’s matched filter from the receiver side into 
the channel. The complexity of the transmitter increases; the 
receiver captures a preequalized CIR which is, in the case of 
large delay spread channels, much simpler to equalize than the 
original one. Research on TR touches a number of different 
research areas. Seen as a preequalizer, TR is known as a 
transmit matched filter [21] or pre-rake [22] which has been 
studied in connection with CDMA but not at very large delay 
spreads, where spatial and temporal focusing comes into effect. 
TR has to be compared to optimal schemes (e.g. [23], [24]). 
Complexity will become a major issue since delay spreads of 
100 or more taps can make many known optimal schemes 
infeasable. The spectral efficiency is also of great interest. 
TR does not use the capacity achieving power allocation [25], 
but in fading channels, an increased bandwidth yields channel 
hardening which could increase outage capacity. Finally, many 
practical issues need to be addressed. At which bandwidth 
does it become advantageous to precode transmitted messages 
in order to shorten the channel, and which precoding method 
should be applied? At which size of a network does it 
become useful to shift complexity from several receivers to 
one centralized transmitter? We believe that the recent interest 
in transmission over very large bandwidths such as provided 
by the standard IEEE 802.15 [26] is a good incentive to rethink 
many of the long known methods that were invented to deal 
with channels of low and moderate delay spread. We propose 
TR as a suitable starting point for these investigations, as it 
clearly points out the potential gains that high-delay spread 
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channels may offer despite their draw backs. 

Here, we shall present an initial study of TR applied to radio 
wireless systems. In Section II, we briefly introduce TR, to¬ 
gether with the formalism used in this paper. In Section III, we 
compare the rate of a TR system to one that uses waterfilling 
and we simulate the cumulative distribution function of the rate 
of such a system in a fading channel. Section IV deals with 
the characterization of time-reversed random fields in wireless 
channels. To this end, we first define suitable metrics in order 
to characterize temporal and spatial focusing. Then we present 
measured results of TR space-time focusing in an indoor Ultra 
Wide Band (UWB) radio channel in the 2-8 GHz band [28]. 
To provide more insight in the space-time focusing results, we 
finally perform simulations for broadband systems (bandwidth 
up to 100 MHz) using two very different channel models at 
2.5 GHz [27]. 

In Section V, we investigate the applicability of TR to 
wireless data transmissions. We examine how TR performs 
compared to optimized linear precoding, and we compute bit 
error rates of both a TR and an optimized system [29]. We 
describe in Section VI how commercially available MIMO 
channel sounders can be converted into sounders that can be 
used for real-world TR-channel measurements, and we finally 
conclude in Section VII. 

II. Time Reversal Formulation 

Consider a transmission between transmit point P and 
receive point Q. The CIR is denoted as h(r,P —> Q). It 
is important to note that the usual symmetry properties imply 
that /i(r, P —+ Q) = h(r, Q —► P). 

The so-called CIR is actually the convolution of the infinite- 
bandwidth physical channel response and the filter impulse 
response of bandwidth B, e.g. a Nyquist filter with given roll¬ 
off factor. 

Depending on the scattering channel and the bandwidth, 
the CIR results in a temporal spreading of the initial pulse. 
Scatterers indeed create multipath mechanisms which, in turn, 
cause echoes to arrive at the receiver with different delays. 
However, the resolvability of the different delays depends 
upon the ratio of the inverse of the bandwidth to the physical 
channel spread (i.e. the interval between successive delays). 
The smaller X/B is relative to the channel delay spread, 
the larger is the number of resolved paths. As an important 
conclusion of the above discussion, it can be expected that the 
B x trms product will play a significant role in TR focusing 
results. 

Using the proposed formalism, the received signal at any 
point R for a pulse emitted from R 0 is h(r,Ro —► R). 
At a particular point T, which we define as the transponder, 
h(r, Ro —* T) is captured. If the transponder sends back the 
time reversed version of the captured signal, i.e. h(—r,R 0 —► 
T), then at any point R, the received signal can be expressed 
as 

s(r, R) = fc(r, T R) * h{-r, R 0 T) (1) 

where * denotes the convolution product. Note that (1) as¬ 
sumes a perfect estimation of h(r, R 0 —> T). In practical 
settings, noise and interference considerations will cause this 


estimation to be biased. It should also be remembered that 
s(r, R) is related to the bandwidth B. 

From now on, we define the point R = R 0 as the focal, the 
target point or simply the receiver. Furthermore, when there is 
no ambiguity about the transmitter position, we also drop the 
symbol T and the arrow in the notation of die CIR. 

Based on (1) and on the symmetry properties mentioned 
above, the signal received at Ro is 

«(r,Ro) = ft(T,Ro)*fc(-r,Ro) (2) 

As a consequence of (2), the time reversal operation causes 
the received signal at Ro to be focused in both time and space 
through constructive interference, meaning that all multipath 
signals add coherently at Ro, and incoherently elsewhere. 

Models usually compute die channel transfer function 
H(uj, Ro) over the system bandwidth (where uj is the angular 
frequency). It is the Fourier transform of /i(r, R 0 ). Relation¬ 
ships (1) and (2) are easily written in the frequency domain, 
since the time reversal operation corresponds to a complex 
conjugation in frequency (denoted by the superscript *): 

Sfa Ro) = H(u>, Ro)JT (u/, Ro) (3) 

According to (3), we observe that the time reversal operation 
is equivalent to a perfect channel matching, i.e., the channel 
acts as the system's matched filter. According to (2), we see 
that s(t, Ro) is the autocorrelation function of the channel, so 
that there is a second interpretation for TR. We may consider 
the CIR as a randomly generated code sequence which is 
given by nature. Then, a TR communication system operates 
similar to a broadcasting CDMA system with a pre-RAKE 
[22]. But the code sequence is not a spreading sequence 
in the sense that it spreads transmission bandwidth without 
increasing information rate for each single user. In TR, the 
fully spread bandwidth can be used for transmission to each 
user. However, since the sequence does not have carefully 
designed correlation properties but is at best random, the 
TR system cannot erase intersignal interference (ISI) and 
multiuser interference completely. 

III. Capacity and Outage in TR Schemes 

We compare in this section the rate of a TR system to one 
that uses waterfilling and simulate the cumulative distribution 
functions of the rates of these systems in a fading channel. 
TR enforces a particular power allocation scheme, P(lj) ~ 
\H(iv)\ 2 , which is in general not capacity achieving. The 
power allocation strategy inherent to TR is to enhance the 
strong components of the channel, and to weaken those which 
are less strong. At low SNR, this strategy is related to water¬ 
filling, so that we can expect a close-to-optimal performance. 
At high SNR, where a constant power allocation scheme is 
optimal, TR is not. We demonstrate this behavior in Fig. 1 
in which the mean rate of a TR communication system that 
operates in a channel with 50 i.i.d. Gaussian distributed taps 
is depicted (dashed line). The rate equals at low SNR the 
capacity of the channel (continuous line), and deviates more 
and more at higher SNR Also for comparison, the capacity 
of the AWGN channel is shown. It is largest in this picture 
since all curves are normalized to a bandwidth of 1 Hz [25]. 
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IV. Time Reversal Channel Characterization 
A. Characterization Metrics 

The space-time received signal s(r,R) is converted into 
performance metrics by considering 

• the energy of s(r, R) at any point R in space, at a given 
time ro, i.e. 

77 £>(R) = |s(r 0 ,R)| 2 (4) 

with ro such that |s(r 0 ,Ro)| = max r {|s(r,Ro)|}. 

• the RMS delay spread of s(r, R) (on a realization basis), 


Ar(R) - 


7 (- 


1 m 


) 2 |s(t, R)| 2 dr 


/|s(r,R)| 2 dr 



Fig. I. The average rate of a TR communication system reaches at low SNR 
the capacity of the frequency selective channel. 


where r m is the average delay defined as 

TVn(R) = f r|s(r, R)| 2 dr/ J |s{r,R)| 2 rfr. 


Cumulative Distribution 



Note that the above definitions assume that the received 
space-time signal is sampled correctly in time so that the 
maximal amplitude can be captured. An alternative definition 
of T)( R) would be to consider the peak energy, irrespective of 
the time delay: 

Vm( R) = max{|s(r, R)| 2 }. (6) 

T 

Both t]d (R) or VM (R) ami Ar(R) are random spatial func¬ 
tions, which can be characterized by their first-order moments. 

1) Spatial Focusing: The spatial focusing around the focal 
point is described by two parameters. 

The asymptotic space-focusing gain (SFG) is given by 

Pd = iFf »7 £>(Ro)A?d(R)- (7) 

|R—Ro |—*>00 


Fig. 2. In channels with a larger number of taps, the outage probability 
decreases. 


In fading channels, outage capacity can be more relevant 
that ergodic capacity. The cumulative distribution function 
(cdf) of the rate of a TR system at 10 dB SNR in channels 
with 10, 50, and 100 i.i.d. Rayleigh taps is shown in Fig, 2. 
The cdf of the TR system is shown as the blue curves, those 
of a scheme that uses water filling in black. The rate of the 
TR system has for a particular number of taps a cdf which 
is very similar to that of the optimal scheme but shifted by a 
fixed amount of less than a quarter bps/Hz to the left. A TR 
system in a channel with 50 taps has at an outage probability 
of 20% about the same rate as an optimally precoded system 
that uses water filling in a channel with 10 taps. The increase 
from 50 to 100 taps does not improve the rate at this outage 
probability much. 


It is the ratio between the energy at Ro to the energy at long 
distance from Ro. A large value of this ratio indicates better 
space focusing. Note that pm is defined similarly with respect 
to 77 m (R), and that pd could be defined at different time 
delays, i.e. other than ro. 

The 3-dB contour of the energy function po (R) or 77 m(R) 
can be considered as the focusing region. It is described by 
the distance in both range and cross-range for which ^(R) 
or 7 /m (R) remains within 3 dB below the energy at Ro. The 
characteristic parameters G a and G x are therefore defined (e.g. 
for rjo such that: 

7?d(Ro + G a u a )/77jo(Ro) = 0.5 (8a) 

//^(Ro 4- G x u*)/7?£>(Ro) = 0.5 (8 b) 

where u a and u* are unit vectors, respectively in the 
range and cross-range directions (i.e. in directions respectively 
parallel and orthogonal to the T — R axis). 
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2) Time Focusing: The time focusing at the focal point is 
described by the RMS delay spread of s(r , R) at R — R 0 , 
denoted as Aro = Ar(Ro). Note that this delay-spread is 
expressed in (5), and accounts for the pulse width. 

Finally, a time focusing gain (TFG) is also suitably defined 
by the relative increase of RMS delay spread At(R) at any 
point R compared to Rq. The parameter is denoted as 


a(R) = 


Ar(R) - Arp 
Ar 0 



The asymptotic TFG is given by limj R _ Ro |_ 00 o( R). A larger 
TFG indicates better temporal focusing in the sense that the 
time compression at the focal point with respect to any point 
away from the focal point becomes larger. 

The performance measures descrived above are random 
variables, depending on the channel realization. Hence, they 
should be characterized by their mean and variance taken over 
all realizations. 


B. Measured UWB Channels 

In this first part, we present CIRs obtained from TR as com¬ 
puted from real data measured with a very large bandwidth. 

1) Experimental Set-Up: Measurements were conducted by 
Intel Corp. at off-peak hours to ensure channel stationarity. 
The environment is an office space (40 m x 60 m) with 
many cubicles. Measurements span the bandwidth 2-8 GHz 
with 3.75 MHz frequency resolution. Antennas are vertically 
polarized. The virtual grid on which the receiver is moved has 
a distance of A m /4 where A m is the wave length of the mid 
frequency of the measurements (5 GHz). The antenna is moved 
to a different location with a precise robotic positioner. At each 
antenna position, the channel is measured over a 6 GHz band 
from 2.0 GHz to 8.0 GHz using a vector network analyser. 
The measurements are corrected to compensate for the system 
components (including cable, gain stages, and antennas). The 
height of the transmit antenna is about 2.5 m and that of the 
receive antenna is lm above the floor. 

The CIR is computed by taking the inverse DFT of mea¬ 
surement data. 

2) Experimental Results: We evaluated r\u for 11 different 
scenarios. Hereafter, we only consider the two extreme cases 
which have the worst focusing and yet they are very similar 
to the outcomes for the other scenarios. The first case is 
an essentially line-of-sight (LOS) scenario with a transmitter- 
receiver separation of 3 m. The second case has a separation 
of 11 m and is a typical non-LOS (NLOS) situation. 

A 3D plot of 77 jvt (R-) of these two cases is shown in Fig. 3. 
In this figure, the square grid spans a region of 7A m x 7A m 
where X m ~ 6 cm. We see that spatial focusing works fine 
in both scenarios. In neither of them, however, the peak is 
isotropic; both peaks have one direction in which they fall off 
faster, and another, in which the decay is slower. The structure 
of the peak carries some information about the geometry of 
the environment, i.e., about directions which show faster and 
those which show slower decorrelation in space. 

The signal power level is at least 10 dB lower at a distance 
of 7A m than its value at die receiver. This demonstrates die 
spatial decorrelation very well. 



(a) LOS scenario. 



(b) NLOS scenario 


Fig. 3. Spatial Focusing tja*(R). One shot spatial field realizations for a) 
the line-of-gight b) the non-Iine-of-sight scenario. 


In Fig. 4 we illustrate the time compression of the received 
pulse. Displayed are the magitude of the CIR for the LOS 
scenario (a), the respective time-reversed, i.e., compressed CIR 
(b), the CIR for the NLOS scenario (c), and the compressed 
one for this case (d). The channel’s power that is not in 
one (central) tap is significantly shortened. Yet, for the LOS- 
case, comparably strong side lobes remain visible, since the 
deterministic component causes correlation. The CIR of the 
NLOS channel is much better compressed in time. Still, 
there is some intersymbol interference in the time compressed 
field and this would require equalizing the received signal. 
However, the time compressed channel has a significantly less 
power in its side lobes and is thus less expensive to equalize. 

For these two sets of data, we computed the temporal and 
the spatial focusing gain which are displayed in Tab. I. The 
TFG is 75% and 30%, respectively, for the LOS and the NLOS 
case. The CIR with file stronger deterministic component has 
a higher TFG. The SFG is 21 and 42 so that the SFG in the 
LOS scenario it is lower. A dominating strong component in 
the channel yields higher TFG, but lower spatial orthogonality. 

To provide a better physical understanding of these observed 
results, we will present in the following a detailed character- 
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TABLE I 

TFG AND SFG COMPUTED FROM EXPERIMENTAL DATA. 


TFG (eq. (9)) 

PM (eq. (6), (7)) 

LOS 0.75 

21 

NLOS 0.31 

42 


or urban canyons. 

Here multipaths are caused by two different mechanisms: 
reflections from the waveguide walls and scattering by the 
medium inhomogeneities. Fig. 5 shows a sketch of the waveg¬ 
uide: two rays are propagating from the source at R 0 to the 
transponder T. One ray reaches T after several reflections 
from the waveguide walls. The other one reaches T by 
multiple scattering on the inhomogeneities. In fact the source 
at Rq produces an entire angular spectrum of rays over the 
bandwidth, and signals reaching T undergo a mixture of these 
two multipathing effects. 



iLl ill ..... 


500 lOOO 

taps 


1500 


(C) NLOS 


.42 0.6 
=! 0.5 
® 0.4 



lOOO 1500 

taps 


2000 


(d) NLOS 


Fig. 5. A sketch of the random waveguide mode! with two rays from the 
source at Ro propagating to the transponder at T each undergoing one of 
the two mechanisms for multipathing: refection from walls and scattering by 
inhomogeneities. 


The waveguide width is L x and the distance along the 
waveguide between Ro and T is L z . We assume that the 
initial transmitted signal at Ro can be given a limited angular 
aperture d. For the sake of simplicity, the mean refractive index 
is chosen as unity (note that changing this value only affects 
the overall phase of the computation). The fluctuation is an 
isotropic, Gaussian correlated random function of space, with 
RMS height ms and correlation length i. 


Fig. 4. Magnitude of channel impulse responses (a and c) and the time 
compressed impulse response (b and d) 

ization of TR focusing in two different wireless channels at 
2.5 GHz: a randomly-filled waveguide and a discrete scattering 
channel. 

C. Random Waveguide 

1) Channel Model: We investigate in this section guided 
wave propagation in two dimensions. The interior of a periodic 
waveguide is a weakly scattering medium. Signals propagating 
through this waveguide experience significant multipathing 
yielding a large delay-spread. This model can seen as an ap¬ 
proximation of propagation through hallways inside buildings 


To compute the random space-time field, we place equi- 
spaced phase screens inside the waveguide separated in dis¬ 
tance by L Simulating propagation along the waveguide in¬ 
volves combining a sequence of steps from one screen to 
the next. Each step involves propagation through an “empty” 
waveguide followed by a random phase correction. 

We now describe a single step of the phase screen method. 
Let z n -i, z n and z n +1 denote locations of three consecutive 
phase screens. We assume that the channel transfer function 
z n —\ —► z n ) from to z n is known. Because the 
waveguide is periodic, we express it as a Fourier series: 


->• z n ) - a m {u>,z n )e 3 ' 2lTmx/Lx . (10) 

m=—oo 
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The channel transfer function for z n to z n + x is then given by 


H(u,z n —► z n +i) = e ^A(x)<5z/ 2 

oo 

J2 T m (u>,6z)a m (u>,z n )ei 2 * mx ' L ’ 




(ii) 


with 6 z = z n+ i — z n and k = 2n/X the wavenumber. The 
forward propagation operator T m (u, Sz) is defined as 


T m (cj, Sz) = exp 



(27 vm/(kL x )) 2 Sz 


( 12 ) 


The random phase correction £l{x) is a path-integral of the 
fluctuation /x(x, z) = 1 — n 2 (x, z) between the screens z n and 
z n +1 with n(x, z) denoting the index of refraction inside the 
waveguide. 

To limit the angular aperture for the initial transmit at Ro, 
we filter the Fourier modes before updating the field to the 
next phase screen. Let zq denote the phase screen containing 
Ro and '&(uj,zq,x) denote the initial transmit signal which 
we express as the Fourier series: 


oo 

y(u,z 0 ,x)= ^ ■4> m {u})e> 2lTrnx/L *. (13) 

m——oo 


From (12) we determine that propagating modes are those for 
which 



(14) 


All others are evanescent. For a propagating mode, the prop¬ 
agation angle is determined from sin$ m = 27rm/(kL x ). The 
largest value of m such that (14) holds gives the largest 
propagation angle $max- To set the angular aperture of the 
initial transmit we restrict the summation in (14) to modes 
m such that < d. 

2 ) Simulation Results: Numerical simulations in the waveg¬ 
uide have been carried out using the following parameters: 

• the central frequency is 2.5 GHz 

• the length L z is 100 m 

• the refractive index fluctuation has RMS height /irms — 
0.05 and correlation length i — 1.0 m 

• the waveguide width L x is either 12.8, 25.6 or 51.2 m 

• the source angle d is either 25°, 45° or 65° 

• the bandwidth B is between 10 and 100 MHz. 


The focus point Ro is centered in the waveguide. We set 
Ro to be the origin for the x — z coordinate system. The 
transponder is located 100 m down-range from Ro. It is an 
aperture spanning 1 m < x < 2 m. Note that the aperture 
center is deliberately slightly skewed from x — 0 as to break 
any symmetry. Hence, any observed effect is by no means 
an artifact caused by symmetry. Note also that our numerical 
results have shown that the space-time focusing is insensitive 
to the location of the transponder array as long as there is 
sufficient multipathing in the channel. 

Figs. 6 and 7 show t?£>(x, zo)/r)i>(xo,z 0 ) and A t(x,z 0 ) 
respectively for a waveguide of width 25.6 m. Averages were 
computed for 500 realizations of the random refractive index. 
Results are shown for bandwidths of 10, 20, 40, 80 and 100 
MHz. 



Fig. 6. Energy of sb{tq,x ) in a waveguide with width L x — 25.6 m and 
L z = 100 m averaged over 500 realizations of the random refractive index. 
The source is limited in angle by A# = 45° and in bandwidth by 100 80 
40, 20 and 10 MHz. 



Fig. 7. RMS delay spread trms averaged over 500 realizations of the random 
refractive index. All other parameters are the same as in Fig. 6. 


In Fig. 6 we observe that quality of spatial focus about x 0 = 
0 is nearly the same over all bandwidths. However, the energy 
away from the refocus location decreases as the bandwidth 
increases. In particular, we observe a 5 dB difference between 
the 10 MHz and 100 MHz cases. 

In Fig. 7 we observe the quality of temporal focus about 
Xq ~ 0 where Ar is the smallest. Because a larger bandwidth 
yields a shorter pulse, Ar is smaller for larger bandwidths. 
Nonetheless, the curves become smoother as the bandwidth 
increases. This phenomenon is due to the onset of statistical 
stability manifested from broad bandwidth signals [34]. 

To examine spatial focusing in greater detail, we show the 
space-focusing gain in Fig. 8. The top plot is for L x = 25.6 
m with different values of $. The bottom plot is for t) = 45° 
with different values of L x . As we have already mentioned, 
adding bandwidth to the system reduces the energy away from 
the refocus location. Fig. 8 further demonstrates this since all 
curves show that SFG increases with bandwidth. The top plot 
in Fig. 8 shows that limiting reduces the SFG. The bottom 
plot shows that there is little difference between the L x ~ 12.8 
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Fig. 8. Asymptotic space-focusing gain (SFG) as a function of bandwidth 
for different source angles (top) with waveguide width L x — 25.6 m and 
different waveguide widths (bottom) with source angle At? — 45°. 


m and 25,6 m cases. However the SFG for the L x = 51.2 
m case is much smaller than the other two. For a fixed d, 
widening the waveguide reduces the diversity of modes that 
propagate in the channel This reduction, in turn, reduces the 
quality of focus in time reversal 
In addition, we show the 3-dB contour in cross-range G x 
in Fig. 9. Similar to Fig. 8, the top plot is for L x = 25.6 
m with different values of t? and the bottom plot is for 
i? = 45° with different values of L x . Here we observe 
that cross-range dependence of the 3-dB contour changes 
very little with respect to bandwidth. Because limiting the 
angular aperture of the source is a spatial low-pass filter 
yielding broader initial transmit sources, G x increases as # 
decreases. For different waveguide widths, the 3-dB contour 
does not change significantly. It decreases slightly as the width 
increases because of added angular diversity. From (14) we 
understand that the number of propagating modes within a 
fixed angular aperture is set by A and L x . With A fixed and L x 
increasing, more propagating modes exists inside the angular 
aperture. For example, with t? = 45°, there are approximately 
75, 150 and 301 propagating modes for L x = 12.8 m, 25.6 m, 
and 51.2 m, respectively. Time reversal methods exploit this 
added spatial diversity to tighten the refocusing. 



bandwidth (MHz) 

Fig. 9. 3-dB contour in cross-range G x as a function of bandwidth for 
different source angles (top) with waveguide width L x = 25.6 m and different 
waveguide widths (bottom) with source angle At? = 45°. 

To examine temporal focusing more closely, we show the 
asymptotic limit of the time-focusing gain in Fig. 10. The 
top plot is for L x — 25.6 m with different values of t? 
and the bottom plot is for i? = 45° with different values of 
L x . All curves decrease monotonically with bandwidth. The 
main reason for this decrease is that the pulse is narrower 
as bandwidth increases. The overall delay spread is reduced 
and the refocus region is not as pronounced. The TFG in¬ 
creases with $ since the additional diversity allows for greater 
temporal focusing in the sense that ratio between Ar at and 
away from xq becomes greater. Similar to the result for the 
SFG, we observe that the TFG for L x ~ 12.8 m and 25.6 m 
are close to each other while the widest waveguide yields a 
smaller TFG. As L x increases effects manifested by reflections 
from waveguide walls decrease and the waveguide domain 
approaches an unbounded one exhibiting a smaller TFG. 

D. Geometry-Based Stochastic Channel 

1) Channel Model: For discrete scattering media, such 
as typical outdoor wireless environments, the channel can 
be represented by a number of effective scatterers randomly 
distributed in space. A ray-based approach can then be used to 
describe die channel as a sum of so-called scattered or reflected 
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Fig. 10. Asymptotic time-focusing gain (TFG) as a function of bandwidth 
for different source angles (top) with waveguide width L x — 25.6 m and 
different waveguide widths (bottom) with source angle = 45°. 


contributions. In the following, the statitical distribution of 
scatterers is based on the well-known one-ring model [30]. 

We first consider that only single scattering occurs. We 
assume that the scatterers are uniformly distributed inside an 
annular region surrounding the target point, as illustrated in 
Fig. 11. This region is specified by a minimal radius p m and 
a maximal radius pM . The angle-spread of the channel (as 
seen from die target point) is fixed by limiting the aperture of 
the annular region to a given portion, i.e. specifying minimum 
and maximum angles, T? m and (see Fig. 11). In this paper, 
we refer to At? = $m — T? m as the scattering angle-spread or 
simply the angle-spread. 

We further simplify the channel description by choosing 
Pm ~ 2 p m and pm « D, with D denoting the distance 
between the transmit and target points. Hence, the RMS delay 
spread of the physical channel (over an infinite bandwidth) 
is well approximated for an omnidirectional distribution of 
scatterers (A# = 27 t) by 


trms « 


2c 


(15) 


with c — 3 x I0 8 m/s is the speed of light. For a limited angle 
spread, we assume that — 7r—Ai?/2 and ~ ?r-bAt?/2. 



Fig. 11. Geometrical representation of the propagation model. 


In that case, the delay-spread can be written as 

y/bpm I At? + sin At? 


trms 


2c V 2 tt (16) 

The channel transfer function is easily estimated by means 
of a ray-tracing approach that yields 


Ro) 

N 

= ^G(u,R 0 ->S i )r i 

dp* s,rfs,T 1 7/2 


1 =1 


N 


= £ r ‘ 


dn^Si + ds t T] 
exp [—jo^dRoS, + <* 8 ,t)] 


G(ui, Si —► T) (17) 


(18) 


1 


[dUoSi + ^SiT ] 7//2 

The quantities in (17) and (18) are 

• u = 2tt/A, the angular frequency (A denoting the 
wavelength); 

• (7(u;, P —► Q) = exp(—ju;dpQ)/(dpQ) 7 / 2 , file spread¬ 
ing function for a transmission from P to Q, with 
7 denoting the path-loss exponent (7 = 2) and <fpQ 
denoting the distance between P and Q; 

• Fj, the scattering coefficient of the 2 th scatterer; in this 
paper, it is considered to be a complex Gaussian variable, 
with given standard deviation; 

• usua l spreading factor for reflec¬ 
tion [^2], ensuring that the transfer function is inversely 
proportional to the total distance, as outlined by (18); 
multiplying this reflection spreading factor by T* yields 
the so-called scattering cross-section. 

Note that (17) and (18) implicitly assume that the CIR 
is filtered by a Nyquist filter with roll-off factor equal to 
zero, hence the filter is a rectangular window in the frequency 
domain. Simulations were done with the following parameters: 

• the central frequency is 2.5 GHz; 

• the transponder is separated from the target point by a 
distance of 10 4 wavelengths, i.e. 1.2 km; 

• when the delay-spread is not used as a variable or 
explicitely assigned a specific value, it is fixed at 0.8/is 
(i.e. p m = 1800A); 

• different angle spreads, ranging from 7r/3 to 27r. 

• the scattering coefficient is assumed to be complex Gaus¬ 
sian distributed with standard deviation of 0.25. 
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2) Omnidirectional Scattering Results: Fig. 12 shows snap* 
shots of tjd(R) (the target point is in the center of the 
figure), for two bandwidths: 50 and 100 MHz, two channel 
delay spreads and omnidirectional angular spreading of the 
scatterers. 


»CM - H - 10 ■ 1 . 61 « - B = 10 UHz 



Rang* ntofcw » MMlangtt 


E{pm) = PM = 2.2 r^B 0 - 52 (20) 


E C®j» ^“ ■ K ° + 111 - *->> 

x |J 0 (2w|R - Ro|/A) | a exp (-|R - Ro|/3.2A) (21) 


with 

a = 2.6exp(-TRMs/l-23) 

(22) 

and 

K D (X l/p D 

(23) 


G a /X = G x /\ = 0.4lT^ s 01 £r°- 02 

(24) 


At 0 = 0.083t^s 25 B _1 + 0.93r^| 

(25) 


£{cr(R)} S 0.37{1 - | J 0 (2x|R - Ro|/A)|^ 

x exp(— |R — Ro[/1.5A)} (26) 


Fig. 12. Typical one-shot realizations of time reversed random fields (the 
energy function is expressed in [dB]). 

The impact of both the bandwidth and the channel delay- 
spread is clearly visible (the scales are kept constant from 
graph to graph). For the smallest B x trms product, some areas 
receive the same level of energy as the target point. For die 
largest B x trms product, the space focusing gain is about 15 
to 20 dB. Also, it seems that the instantaneous energy function 
is oscillating as a function of the distance to the focal point 

Fig. 12 shows that the space-time focusing is improved 
when both the bandwidth and the delay-spread are increased. 
Figs. 13 and 14 show £{77/5(R)/7 Jd(Ro)} and E{At(R)} as 
a function of die system bandwidth. The channel delay spread 
is set to 0.8 fis for these simulations. 



10D 3 to CTi'«*rtfh 


Fig. 13. Simulated energy function E{t)d(JI)/t)d(IIo)} as a function of 
range and bandwidth (trms = 0.8/xs). 

The simulated dependencies have been fitted by the follow¬ 
ing empirical formulas: 

E{pd] =Pd = 5.9r^f]gB 035 (19) 


with j3 — 1.3exp(—Here, the bandwidth B is 
expressed in MHz, and die delay-spread trms is expressed 
in /is. 

The alternative energy ratio /t)m (Ro)} can be 

expressed in die same fashion as in (21). 

• In (19), (24) and (25), exponents for B and trms arc 
similar when both variables appear multiplied by each 
other. This observation implies that space-time focusing 
in general is improved by increasing the product B x 
trms* However, that assertion should be nuanced, since 
some terms in (25) and (26), as well as a and {3 are only 
related to die channel delay-spread. 

• When B is sufficiently large (> 10), die space-focusing 
parameters depend upon B x trms- This impact is sig¬ 
nificant regarding the asymptotic SFG. However, the 
dependence of the 3-dB contour width towards B x trms 
is much weaker. Both G a and G x are therefore mostly 
inversely proportional to the carrier frequency only. Note 
that for smaller bandwidths (for which eqs. (19) and (20) 
are not valid anymore), the average asymptotic SFG tends 
to unity. 

• The average delay spread at the focal point Aro results 
from die additive combination of two terms: one roughly 
related to B x trms, and the other one to trms only. For 
large B x trms, the variation of Ato with B becomes 
weak, so Atq only depends upon the channel delay- 
spread. This is clearly visible on Fig. 14. At any other 
point R outside the focusing region, the average delay 
spread 1 ?{At(R)} is an oscillating increasing exponen¬ 
tial function of the distance from the focal point. It 
saturates at long distances. 

• At any point R, the time-focusing gain £’{(j(R)} is 
only related to the channel delay-spread. Furthermore, 
the asymptotic TFG (i.e. at long distance) is a constant 
value (~ 0.37) independent from B and trms, at least for 
sufficiently large B x trms (> 10). 
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Fig. 14. Simulated RMS delay spread £{Ar(R)} as a function of range 
and bandwidth (trms — 0.8/xs). 


3) Impact of Reduced Angular Spread: When the angle 
spread is decreased from 2ix to smaller angles, the spatial 
focusing is degraded as expected. 

» For Ai? > 7r/3, the asymptotic SFG po decreases as At? 
decreases, and is found to vary as (Ai?/27r) 0 ' 6 . 

• The 3-dB contours G a and G x are not identical any¬ 
more, and G a — 4 G x for At? < 37r/4. They are 
decreased respectively by Ai? -1-3 and At?” 12 relatively 
from their previous values (i.e. when At? was equal to 
27r). For At? > 37 t/4 , die ratio (5 a (At?)/G a (At? = 2n) 
continues to decrease following a At? -1 - 3 law while 
(5 z (At?)/<5 a; (Ai? = 27r) remains constant. Hence, G a = 
G x for At? = 27t. 

New empirical egressions can therefore be derived easily 
that take into account that the channel delay spread trms is 
also modified when reducing the angle spread, as highlighted 
by (15) and (16): 

PD = 5.9r^|B°- 35 (At?/27r)°- 6 I /- a2 
with v — (At? + sin At?)/(27r), 

Ga! A = 0.41r^ s 01 B' 0 02 (Ai?/27r) _1 - 3 i/ 0005 
« 0.41(Ai?/27r)~ 1,3 

f 0.13rjrwo 01 B~ 0 02 (At?/27r)~ 1 * 2 t/ 0 005 , 

G x /X = \ tt/ 3 < At? < 3 tt/4, 

°- 41r RMs 0l5 ”° ° 2l/0 ' 005 ’ > 3?r / 4 

0.13(At?/2^r)- 1 - 2 , ?r/3 < At? < 3 tt/ 4, 
0.41(At?/27r)^ 1 ’ 3 , At? > 37r/4. 

4) Benefits from Multiple Antennas: The impact of using 
multiple antennas at the transmit point is analyzed in the 
following. At the transmit point, the terminal now consists 
of Mt antennas. Each antenna has a specific location around 
T, denoted as T u . Hence, the focused signal is written as 

s(r,Ro) = 

i Mr 

~ EMr,Ro-T u )*h(-r,Ro-T u ) (30) 

f-f. 



(27) 

(28) 

(29) 


where the scaling factor \jyjMr ensures that the total trans¬ 
mitted power is kept constant, irrespective of the number of 
antennas. 

The benefit from multiple antennas is proportional to Mt 
as far as the SFG is concerned, in agreement to well-known 
diversity results [31]. Indeed, considering a constant transmit¬ 
ted power, the peak energy at the focal point will grow as Mt, 
but will remain the same anywhere else outside the focusing 
region. 

5) Multiple-Bounce Model: Single-bounce models are not 
sensitive to the order of magnitude of T and 7 since they are 
only responsible for a proportionality factor. This is not the 
case for multiple-bounce models in which the average loss 
per interaction, related to the scattering cross-section and the 
path-loss, is a critical parameter. 

The model can be easily extended to consider multiple 
bounce interactions. In such case, the transmitted waves are 
scattered from one to a certain number of times. The calcu¬ 
lation of the transmission channel can be written as a sum of 
N scattered contributions: 

N 

tf(u/, Ro) = J'CttnGfaSi —> T) (31) 

2=1 

where 

• G(uj, S 1 —> T) is the spreading function for a transmis¬ 
sion from the I th scatterer to the target point T, 

• 171 is the scattering cross-section of the I th scatterer, 
proportional to the scattering coefficient Ti [33], 

• Ci is die I th element of vector C, which is the solution 
of a JV-order linear system given by: 

C = A -1 • Co (32) 


with 

• Cqi = G(u;,Ro —► Si) is the channel between the target 
point and the location of each scatterer, 

* A is the system matrix accounting for the transmission 
between each pair of scatterers: 


Aki = { 


V 


r k G(u>,s k 

1 


Si) 

k = L 



As already mentioned, this model is highly sensitive to the 
average loss per interaction, i.e. to the product \r) k G{uj^ S k 
S|)|. On one hand, if the latter is very small, then the 
delay-spread of the channel tends to be infinite, which is 
unrealistic. On the other hand, if the loss per interaction is 
very large, the impact of high-order scattered contributions 
becomes small, and die received field is very similar to the 
first-order field. A reasonable assumption is that contributions 
up to the third order can be significant. The average loss per 
interaction should then be chosen accordingly. For an average 
loss per additional reflection of 12 dB, the impact of multiple 
reflections on the focusing is negligible. If die average loss 
per interaction is decreased by 4 dB, the impact, though not 
negligible, remains small. This is mostly because the large 
number of scatterers enables the time-reversed field to be 
highly focused with single-bounce interactions. 


6) Comparison with Random Waveguide Model: In both 
models, the space-focusing gain increases with bandwidth. In 
the geometry-based model, it is also shown that the SFG 
increases similarly with the channel delay-spread. In the 
waveguide model, simulation results illustrate that the SFG 
also increases with the channel delay-spread. The increase for 
the waveguide is related to the source angle and to the inverse 
of the waveguide width. Time-focusing gain examined using 
the 3-dB contour shows a weak dependence on bandwidth 
in both scenarios. Similarly, the dependence on the channel 
delay-spread is also weak, but the 3-dB contour decreases as 
the angle-spread increases. This is clear in the geometry-based 
model. In the waveguide model, we have observed that the 
dependence toward the waveguide width is weak, while there 
is a noticeable dependence on the angle-spread. 

Some differences between these two models have been 
observed. The most pronounced one is that the SFG appears 
to be independent of the bandwidth in the case of the discrete 
scattering model, while a significant dependence is found in 
the waveguide. Also, the absolute value of the SFG is much 
larger in the waveguide scenario. To better understand this 
difference, we need to refer to the specifics of the models. 
In the waveguide case, two mechanisms cause delay-spread: 
large-scale multipath results from the coherent reflections by 
the waveguide walls, small-scale multipath is caused by the 
non-homogeneities of the refractive index. By contrast, the 
geometry-based propagation only consists of non-coherent 
reflections, similar to the waveguide small-scale multipathing. 
Focusing in time reversal is enhanced for bounded domains in 
that it is more robust This accounts for the difference seen in 
these two models. As the waveguide width approaches infinity, 
it begins to agree with an unbounded medium. 

7) Space-time focusing in time-reversed transmissions: So 
far, we have assumed that a single time-reversed pulse is sent 
from the transmitter to the target point. However, in Section V, 
we will deal with data transmissions in TR schemes. As a 
preliminary study, we analyze here how space-time focusing 
is affected when sending a time-reversed data signal. The latter 
consists of the convolution of the time-reversed pulse with a 
train of information bits. 

Fig. 15 shows tiie SFG as a function of the relative data 
rate (i.e. the data rate relative to the pulse bandwidth) for two 
different Nyquist roll-off factors. The bandwidth is equal to 40 
MHz and the channel RMS delay-spread to 0.8/xs. For these 
values the single-pulse SFG pu ~ 13 dB. It is clear from 
Fig. 15 that for the full rate the space focusing advantage is 
totally destroyed. The same conclusion applies regarding the 
TFG, as depicted in Fig. 16. In other words, for a single user, 
the time compression effect resulting from TR is significantly 
reduced if the rate is not seriously decreased. 

In fact. Fig. 15 and 16 show that space and time focusing 
are slowly restored as the data rate decreases. As an example, 
the SFG is about 5 to 6 dB for a rate of B/10, while the TFG 
is around 0.22 for that same rate. For higher roll-off factor, the 
degradation is less than for the zero roll-off case. This is not 
surprising since higher roll-off factors yield lower side lobes. 

Therefore, there is a fundamental trade-off between the data 
rate and the space-time focusing effect that can be achieved 



Fig. 15. Spatial focusing gain as a function of data rate and roll-off factor. 



Fig. 16. Time focusing gain as a function of data rate and roll-off factor. 

in a wireless channel for a given bandwidth. This trade-off is 
dictated by the coherence bandwidth of the channel transfer 
function and is analyzed in more detail in the next section. 

V. Communications over Large Delay Spread 
Channels with Very Simple Receivers 

In this section, we investigate the performance of a most 
simple TR communication system. We assume that the channel 
shortening achieved by TR is all that is done to combat 
remaining ISI in the channel and investigate an uncoded 
MISO system that employs a simple one-tap receiver without 
any further equalization. We study the performance of this 
TR system in the limit of very low and of infinite delay 
spread. We compare this performance to that of an optimized 
preequalizer, and examine the impact of the number Mt of 
transmit antennas, and of a strategy called rate back-off that 
can combat remaining ISI in the channel in a very simple way. 
The goal of this study is two-fold. We obtain a lower bound on 
the performance of a TR communication system in channels 
with very high delay spread, since any other system that does 
employ coding or equalization at the receiver will perform 
better. Note that the upper performance bound was presented 
in Section III, since the computation of the system’s rate in 
that section assumed both infinitely complex transmitter and 
receiver and perfect coding. Secondly, the proposed setup can 
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Fig. 17. A block digram of a MISO system with a preequalizer at the 
transmitter and a single tap receiver. 


be of interest by itself for systems that operate with very large 
banwidths, not too high data rates, and that rely on extremely 
simple receiver architectures. 

Fig. 17 depicts a block diagram of the system. A symbol 
x is upsampled by a factor of D , put on the Mt branches 
to the transmit antennas, filtered with a function g z which can 
possibly be a TR prefilter, and sent over the channels hi. After 
receive noise is added, the signal can be enhanced by a gain 
3, and it is downsampled before detection. 

For all plots shown in this section, we use a Rayleigh fading 
channel with L—100 i.i.d. taps. 

A. Optimal Prefiltering 

First we derive the optimal prefilter for the system in 
Fig. 17. We describe the input-output relation of the system 
in matrix form as 

y[k] = /3xfcHg 4 - 0n[k] (34) 

where x* = [x\X 2 - • • x 2 l-i] is a vector of input symbols, 
y[k] the output symbol and n[k] the noise. The discrete time 
channel matrix H is block toeplitz; its block-columns are 
comprised of the shifted channel vectors [hi [1 ],..., h{\L [], and 
in each block columns that belong to adjacent antennas are 
next to each other. The vector g is given by 

g=[g[0f ... g[£ - 1]T (35) 

and 

g[f] = [gi{l] ... 9mM T l = 0,...,L — l. (36) 

The product Hg is die convolution of the prefilter and the 
channel; T denotes the transpose. 

We give here the optimal g and 3 for the case of D — 1 
(no upsampling). In order to simplify receiver processing, 
the purpose of the set of transmit filters g is to preequalize 
the signal. This means that at the receiver, as much power 
as possible is focused in the central peak, and that ISI is 
minimized, subject to the power constraint || g ||= 1 where 
|| . || denotes the norm. Such a problem can be stated as an 
optimization problem, whose solution is 

g = i (h"H + £l) H"e A (37) 

/J 2 = e^H fH H H + 

A is the delay of the equalizer, and the unit vector with 
1 in entry A, and P Tl = Z?[x£x fc ]. Eq. (37) defines a filter 




Fig. 18. The focused CIR for the TR- and the optimal prefi Iter. 


with an MMSE-type solution with a finite number of taps. We 
notice that the optimized value for /3 bascially ensures that the 
power constraint of the filter is valid. In the transmit-receiver 
chain of Fig. 17, 0 has no effect on the SNR at the receiver; it 
ensures, however, that the transmitter can take into account the 
SNR to optimize the preequalizer. We note that even though 
there exists a striking similarity to an MMSE receiver, we 
could not establish the convexity of the optimization problem. 
Hence, our solution (37) may actually only be a local optimum. 
The similarity of the solution to the MMSE receive equalizer 
suggests that the found optimum shows good performance; but 
better solutions may exist. 

The TR-prefilter is a matched prefilter; it follows from the 
optimal filter (37) for low SNR, 


g TR 


lim 

SNR- 


0 


g — h* 


s.t. h |— 1 , 


(39) 


or 

with 7 = fj2i=i Ya= a IM^ll 2 ) » >- e -> the fil ter equals the 
time-reversed CIR of the channel, and the receiver receives 
the autocorrelation function of the channel. The power in its 
central peak is maximized, but the ISI suppression is worse. 

A picture of the CIRs that result from precoding with the 
optimal and the TR strategy is shown in Fig. 18. The SNR 
is 10 dB. The optimal prefilter suppresses the ISI better than 
the TR prefilter; close inspection reveals that the power in 
the peak of the CIR is slightly higher in case of TR. The 
improved ISI suppression in case of file optimal prefilter comes 
at computational cost The optimal transmit prefilter requires 
a matrix inversion which can be difficult if the number of 
channel taps is too high or the channel is not sufficiently 
stationary. 

B . Rate Back-off 

Fig. 18 suggests that the application of TR or even of an 
optimized prefilter with a single antenna leaves considerable 
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residual ISI in the channel One method to combat ISI is to 
transmit at lower rate; if the transmission rate is by a factor 
of D lower than the bandwidth of the channel, the ISI is 
on average also reduced by a factor of D, but the useful 
power, concentrated in the peak of the CIR, remains constant. 
Such a strategy is not favorable if high spectral efficiency 
is required. But for systems which operate at very large 
bandwidths, such as ultra-wideband systems, it can be very 
attractive to trade abundantly available transmission bandwidth 
against complexity. 

Rate back-off by a factor of D is equivalent to upsampling 
the signal before transmission, and downsampling it by D 
before detection. The input-output relation of the pre-equalized 
channel with rate back-off factor D , D = 1,..., L, can be 
derived as 

=/3]T * 9 i)[Dl]x[k - 1} + pn{k\ (40) 

2=1 l 


where p — ■ 

We furthermore define the SNR pmfb that would be 
obtained if there was no precoding but the receiver had a 
matched filter and used full channel knowledge as 

Pr.ElZfj, Ei Mill 2 ! 


PMFB ~ 


a J 


(46) 


2) Application to TR Systems: The simplicity of the TR 
prefilter allows to compute the ratio k of the useful power 
against the ISI for an exponentially decaying channel with 
delay spread cr T . One obtains in the limit of large L/o r [29] 

KTR = (l-exp-^-j 

|M|(1 +exp-i) + M t ( 1 - exp —i) 


x - 


2M T exp-£-(l-exp-^) 


(47) 


where x[k] is the input symbol and [k] the output symbol 
received with rate back-off. We see that the rate back-off yields 
an equivalent channel which is sampled only at each Dth 
instant. 


C. Focusing Ratio and Effective SNR 

1) Definitions: We introduce some performance measures 

that we will later use for system characterization. A very 

intuitive performance measure is the focusing ratio k which 

we define as , 

E[Po] 


K = 


E[QW] 


(41) 


where E[Pq] is the expected power in the central peak of 
the preequalized CIR, and E[Q^\ is that of the ISI with rate 
back-off D, i.e., the power in the tails of the preequalized CIR. 
The better the focusing, the larger is k. Unlike the measure 
introduced in (9), the measure here takes into account the 
entire, temporally spread power of the received signal, and 
not just the ratio of the energies of the maximum and of one 
distant tap. 

We can express the power of the received, possibly down- 
sampled, signal as 

£[|y [D1 [fc]| 2 ] = P Tx E[P 0 \+P Tx E[Q^] + a 2 (42) 


where we used i?[a;[T[x[m]] = PTx^—tn], with Pt x denoting 
the transmitter power, and E[n[l]n*[m}} = a 2 S[l - m\ with 
a 2 the noise power and S[l — m] the Kronecker Delta. The 
constants Pq and are given as [29] 


P = 




hi\l}gi[m]d{Dk — m — l\ 


(43) 


.(44) 


We define the effective SNR p e j / as 


PefLD 


PtxE[Pq] 
PtxE[QW} + 0 2 


E[Pq] 

E[Q [ D ]] +1 



Using (47) together with (45) and (46) we arrive at an analyt¬ 
ical expression for the effective SNR of a TR communication 
system with Mr transmit antennas and a rate back-off of D, 

MSP = PMFB -T— J m IE ■ (4g ) 

' «TJ? 

We can derive from this the following limiting cases: 


Mf,D 

Peff,D 


Pmfb for ar T —♦ 0 

PmfbDMt r 
-—— for a. 

Pmfb + DMt 


> oo 


(49) 

(50) 


Eq. (50) can be split into a low-SNR and a high-SNR regime 
as 


Peff,D *= PMFB for Pmfb low, o T large (51) 
Peff.D = DMt for pmfb, v t laige. (52) 

Interpretation of these equations yields the following conclu¬ 
sions: 

• Eq. (49) tells that a TR system is lossless if the delay 
spread of the channel is very small Preequalization of the 
channel does hardly change the performance if the delay 
spread is very small We recall from Fig. 2, however, 
that the system in this case is subject to fading which 
will reduce the rate. 

• Eq. (51) shows that TR is optimal at low SNR; this 
behavior is expected for a transmit matched filter which 
maximizes the power in the peak of the signal 

• Eq. (52) demonstrates that in an environment with a large 
delay spread, the performance of a TR system at high 
input SNR saturates at a value independent from the delay 
spread of the channel The performance increases with 
3 dB per added transmit antenna, and with another 3 dB 
for each rate back-off by a factor 2 from the maximum 
available bandwidth. 

The last property is remarkable, because it holds for an 
uncoded system with no receive equalizer, a fairly simple 
transmit equalization strategy, and arbitrary delay spread, i.e., 
arbitrary bandwidth. In a 1 GHz time-invariant channel, a sys¬ 
tem with a powerful two-antenna transmitter can theoretically 




14 




Fig. 19. The effective SNR increases with the number of transmit antennas. Fig. 20. Rate back-off improves the effective SNR. 


achieve a transmission rate of 125 Mbps at 12 dB receive 
SNR with a single one-tap receiver and independent of the 
delay spread of the channel. 


D. Performance Evaluation 

We study the performance of the optimal and the TR 
prefilters in dependence on the number of transmit antennas 
and the rate back-off, and examine furthermore the variance 
of the received power as a function of bandwidth. 

In Fig. 19 we plot die effective SNR of the two filters 
as a function of pmfb ranging from -10 to 20 dB; in this 
picture, an optimal system would be characterized by a line 
with slope 1. The effective SNR is computed using 100 CIRs. 
The curves for die MMSE filter are dashed, those for the TR 
filter continuous. The three different curves for each filter are 
obtained for Mt = 1,2, and 4. The results for the TR filter 
match the prediction of (47). For large pmfb* the effective 
SNR becomes constant. It is 0 dB for Mt = 1, and increases 
with 3 dB for each additional transmit atenna. At low SNR, 
the curves increase linearly with pmfb * hi this range, the 
TR filter performance equals that of the MMSE filter. The 
effective SNR for the SISO system is quite bad; however, the 
shown performance is achieved with an extremely simple, one- 
tap receiver, and in a channel with very long delay spread. 
The optimized prefilter shows much better performance. For 
Mt = 1, the effective SNR saturates at a value which 
appears to depend on the delay spread of the channel and its 
correlation between taps. For more than one transmit antenna, 
the optimized prefilter can completely eliminate ISI, so that 
p e ff grows without bound. The slope of the curve depends on 
the number of transmit antennas. We believe it is interesting 
to study analytical expressions for the behavior of p e jf for 
the optimal prefilter in the limit of very large bandwidth. 

The system can operate with rate back-off in order to 
increase the performance at the cost of transmission rate. The 
effective SNR for a single-antenna system with rate back-off 
of D = 1,2,5, and 25 is displayed in Fig. 20. It increases with 
decreased rate. Also, the saturation of the effective SNR sets in 


much later; the (almost) linear increase of p e j / reaches for the 
TR system with D = 5 up to pmfb = 10 dB, and for one with 
D — 25 up to 15 dB. In this linear range, the performance 
of the TR prefilter is, in terms of the effective SNR, partly 
higher than that of the MMSE prefilter, since the latter is 
optimized for ISI suppression at D = 1 and has less power in 
its central peak. If the MMSE prefilter had been optimized for 
transmission with rate back-off, it would be at least as good 
as the TR prefilter. The design of such a prefilter could be 
interesting for preequalized wideband data transmission. When 
spectral efficiency has not highest priority, transmission rate 
back-off (while keeping the transmission bandwidth constant) 
is very efficient in order to improve system performance. 

We characterize the fading of the channel, which we com¬ 
pute as the normalized variance of the instantaneous effective 
SNR, 


a SNR ~ Var 


PtxPo 


/ E[- 


PtxPq 


PtxQM + o 2 ' l P T *QW + <r 2 




(53) 


We compute the variance for both filters for p = 10 dB and 
Mt — 1 in dependence on the number of taps of the channel. 
If the channel has only a few taps the peak of the preequalized 
CIR fades with higher probability than at high bandwidth with 
a large number of taps. The result is shown in Fig. 21 for 
the number of taps in the range between about 1 and 100. 
The shown curves are based on an average over 100 different 
realizations of the channel. We see that the variance decreases 
quickly with the number of taps; at only 20 taps it is already 
very low. 


E. Bit-Error-Rates of the System 

We plot the bit-error-rates (BERs) of the system, for rate 
back-offs of D — 1,2,5, and 25, i.e. the same as in Fig. 20. 
In Fig. 22 we show the BER of the TR (blue, cont.) and 
the MMSE (black, dashed) prefilter against SNR. Simulations 
were obtained from sending 5000 bits over 5000 different 
channels with L = 100. For D — 1 and also D — 2, BER of 
the TR prefilter floors quickly. For large rate back-off, BER is 
due to the increased effective SNR in the same range as that 






























15 


* 



Fig. 21. The variance of the received signal decreases with increasing 
bandwidth. 



Fig. 22. The BER of the TR and the MMSE prefi Iter. 

of the MMSE prefilter for lower rate back-off. The TR system 
performs well for D > 5; the MMSE prefilter performs well 
for D > 2. For large D, both systems perform in general 
similarly well, since the MMSE filter was not optimizted for 
rate back-off. A worthwhile study would be the performance 
of a system equipped with a short MMSE equalizer at the 
receiver. The performance would potentially be much better 
and the rate back-off could be reduced even though only a 
simple receive equalizer was used. 

VI. TR Channel Sounding 

The value of TR for practical systems can only be assessed 
through experiments. In this section we describe how a channel 
sounder can be designed that can perform real-world TR 
experiments. This design does not require new technological 
concepts; we rather want to show how a conventional wide¬ 
band MIMO sounder can be converted into a sounder for TR 
experiments. A more detailed version of the following can be 
found in [35]. 


A. Architecture 

In TR experiments, one wants to demonstrate temporal and 
spatial focusing achieved in a wide band channel. In order to 
demonstrate temporal focusing, high resolution of the received 
CIR is required. To reveal spatial focusing, more than one 
receive antenna is necessary; one receive antenna is the one 
that receives the time compressed signal. At any other antenna, 
a noise-like signal will be visible. A wideband SIMO system 
is hence the minimum requirement to demonstrate TR. The 
key difference between between MIMO channel sounding and 
a TR experiment is that the latter requires CSI. The CSI 
must be updated in intervals the length of which depend on 
the Doppler spread of the channel. Since even comparably 
stationary channels as those for fixed wireless access can 
have a Doppler spread of more than 1 Hz [36], the CSI 
must be updated much faster than once per second and slow 
methods, like a frequency domain sweep measurement, are not 
useful. We will consider time- or frequency domain correlation 
methods, since they are both fast and can be used for wideband 
channel sounding. 

In common correlation-based wideband channel sounding 
experiments, a transmitter transmits known sequences to a 
receiver. Transmitter and receiver are synchronized, and the 
receiver uses the received sequences in order to estimate 
the channel. In a TR experiment, the transmitter knows the 
channel, and convolves the transmit sequence with the time 
reversed CIR before transmission. The receiver uses again 
the received sequences to estimate the time-reversed channel, 
which will now exhibit spatial and temporal focusing. In 
order to convert a MIMO channel sounder into a TR channel 
sounder, CSI must be made available to the transmitter, and 
the transmitter must be able to convolve the CIR with the 
transmit sequence. Furthermore, the sounder must still be able 
to estimate the channel before CSI is fed into the transmitter. 

A high-level sketch of a possible design is shown in Fig. 23. 
A MIMO channel sounder, which consists of an RF part for 
transmitter and receiver (RFTx and RFRx), a signal generator 
at the transmitter (SGTx) and a data storage unit at the 
recevier (DSRx) is equipped with a TR-signal processing 
unit (TRSPU). The main task of the TRSPU is to drive the 
estimiation of the CSI, and to transfer the CSI via the feedback 
loop to the transmitter, where the CIRs are convolved with 
the transmit sequence and transmitted over the channel. The 
details of the TRSPU and the feedback loop are explained in 
the next subsection. 

B. TRSPU and Feedback Loop 

A TR channel sounder operates in three steps, which are 
controlled by die TRSPU and shown in Fig. 24. In a first step, 
the channel is estimated between a given transmit antenna, 
say Tx o, and a given receive antenna, Rxq. This can be done 
by sending the known sequence T(uj), here expressed in fre¬ 
quency domain, n\ times over the channel. Once the channel is 
estimated with some predefined accuracy, the estimated CSI 
1Tg,o is sent via the feedback loop to the transmitter. In a 
third step, the known transmit sequences T(w) are multiplied 
with the estimated transfer function and transmitted into the 
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Fig. 23. A TR channel sounder can be designed using a MIMO channel 
sounder with a feedback loop for CSI. 



Fig. 24. The operational steps of a TR channel sounder. 

channel. In this third step, the sounder switches between all 
its receive antennas. At Rx o, it will receive a time reversed 
CIR; at all other antennas, it receives a noise-like signal. This 
time-reversed sequence can be sent ri 2 times, until the CSI 
is outdated and the sounder has to start again with channel 
estimation. 

Particular care is required for the choice of the feedback 
loop. The transferral time of die CSI must be much shorter 
than the coherence time of the channel. The data rate of the 
link becomes hence a critical parameter. A single complex CIR 
with about 1000 taps and 8 bit resolution has a size of about 
16 kbit. Assuming a channel with Doppler frequency of 1 Hz, 
a realistic transmission time should not be longer than 0.2 s; 
consequently, the feedback loop needs a data rate of at least 
80 kbit/s per CIR that is transferred. 

The following options can be considered. 

1) Microwave Link: A microwave link is advantageous in 
terms of data rate, which is in the order of Mbit/s, and in terms 
of flexibility. Microwave links can be established whereever 
the experiment is to be performed. However, TR experiments 
are likely to be conducted in environments with large delay 
spreads. The TR sounder will operate at maximum allowed 
EIRP. This limit holds for the microwave link, as well. In 


areas where the SNR of the TR experiment is low, that of 
the microwave feedback link is low, too. The link quality can 
be improved by using directional antennas mounted on high 
poles, but this requires difficult adjustments and large efforts 
in the preparation of the experiment. 

2) Cellular Systems: Currently available cellular systems 
such as GPRS offer great flexibility without the need for 
heavy experimental equipment. Their main drawback is their 
comparably low data rate (order of 100 kbit/s for GPRS) and, 
in particular, the unreliable transmission delay. Since the data 
rate offered does not have too much overhead compared to the 
data rate required, the application of current cellular standards 
for the feedback loop appears risky in particular if the Doppler 
spread of the channel may at times exceed 1 Hz. 

3) Ethernet: The ethemet offers tremendous data rates; its 
main drawback is the inflexibility of its access. However, most 
TR experiments will be conducted in urban areas, where access 
to the ethemet can be expected to be in range. A wireless 
network system such as one based on the IEEE 802.11b 
standard still offers data rates in the order of Mbits and can be 
employed to bridge the gap between the experimental device 
to the next ethemet connector. 

The final decision about the best link depends on the 
availability of particular services, on local network traffic, and 
the Doppler spread of the channel. 

C Further Critical Specifications and Sources of Error 

1) Range and Power: TR works best in environments with 
large delay spreads, i.e., large distances. Consequently, the 
demands on the link budget are high. Critical parameters 
which determine range are here pathloss exponent, maximum 
peak and average transmit power determined by hardware and 
regulations, bandwidth and carrier frequency, antenna gains, 
correlation gains, the sounder’s sensitivity and target SNR. 

2) Phase Noise: TR requires exact phase information at 
the transmitter. Inaccuracies in the phases occur since a) the 
channel is not perfectly time variant, b) the LO signals derived 
by PLLs from die (usually) 10 MHz reference have phase 
noise. The time variance can be overvome if the data rate of 
the feedback loop is large enough. 

Phase noise destroys the coherent superposition of signals 
transmitted from multiple antennas. Since it enters both at the 
transmitter and the receiver, it will reduce the received power 
by a factor of cos(a<£i + a where a 0i, i — 1,2 are 
random variables that describe the phase fluctuations of the 
transmitter’s and the receiver’s oscillators. For phase errors 
less than 15°, this yields a performance loss of about 10 % in 
SNR, which is tolerable. 

VII. Conclusions and Outlook on Future Work 

In this paper we introduced die concept and goals of time 
reversal (TR) techniques in the context of broadband wireless 
communications. TR by itself is a very simple preequalization 
scheme which moves the system’s matched filter into die 
radio channel. When applied to wideband channels with very 
large delay spreads, this preequalizer very effectively captures 
the particular advantages of this channel, i.e. temporal and 
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spatial focusing as well as channel hardening. TR allows to 
shift complexity from transmitter to receiver and, thanks to 
its comparably low computational complexity, is likely to be 
a feasible scheme even if original channel impulse responses 
are very long. 

We showed that the average rate of a TR scheme at low 
SNR is very close to the ergodic capacity of the frequency 
selective channel. In fading channels, the outage capacity can 
be improved if one transmits at higher bandwidth and exploits 
the channel hardening provided by TR. 

Temporal and spatial focusing were investigated with three 
examples based on mesaured indoor UWB data, and on models 
of discrete and of continuous scattering. We illustrated how 
several proposed characterization metrics vary as a function of 
system and channel parameters, such as the bandwidth, delay 
and angle spreads, and number of antennas. 

When addressing communication theoretic issues, only the 
most simple communication scheme with a single tap receiver 
has been considered. We saw that such a system can work in 
this harsh channel when the transmission bandwidth is sacri- 
fied for temporal/spatial focusing. More optimal transmission 
methods can achieve much better performance at the cost of 
more complexity at the transmitter. 

Finally, we outlined how the TR channel can be measured 
with a system which is basically an upgraded commercially 
available MIMO channel sounder. 

This overview touches many issues about transmission 
over channels with large delay spread for which TR is one 
potentially useful transmission scheme. Many questions have 
not been addressed. TR is optimal at very low SNR for very 
simple receivers. A natural question is to ask what happens if 
the receiver has more than just one tap, but still far not enough 
to equalize the channel on its own. It is not obvious that this 
problem results in a MMSE receiver, which could be used in 
combination with the TR precoder. 

No space in this paper has been devoted to the analysis of a 
communication system that makes use of the spatial focusing 
properties of the proposed preequalization schemes. This topic 
is inherently connected to multiuser communications. The rate 
region of such a system requires the study of fading, frequency 
selective broadcast channels in die limit of large bandwidths. 
The derivation of optimal transmit and receive schemes would 
continue the work presented in this paper, or related work 
about joint transmit/receive optimization for multiuser commu¬ 
nications but in the case of large delay spread. For TR, there 
appears to be a natural relation to the interference channel. In 
TR, the channel impulse response is temporally and spatially 
compressed solely with the knowledge of the channel of one 
individual link; no knowledge about the link to any other user 
or network is required. Hence, another potential application 
of TR is wideband networking with strong access points, 
architecturally very simple users, and a very low degree of 
cooperation. 

A very critical aspect that has only briefly been discussed 
in Section VI is that of channel estimation. We believe that 
TR or any other preequalized system can be applied in TDD 
systems where the channel state information is estimated at the 
base station in the uplink, and then used in the downlink. That 


way, only the powerful base station needs to fully estimate the 
channel, and the handsets only estimate as many coefficients 
as the short length of their receive equalizers requires. In the 
uplink, no preequalization is performed, and the matched or 
any other receive filter is not in the channel, but in the receiver 
itself. 

We believe that the study of communications over large 
bandwith / high delay spread channels is very fruitful. These 
channels offer potentially large rate, high diversity, and spatial 
and temporal focusing. TR systems that operate here provide 
a low-complexity solution for data transmission; their inves¬ 
tigation may be a suitable starting point for more research in 
this area. 
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Abstract 

Using broadband radio wireless measurements in an indoor 
environment we demonstrate the remarkable space-time focus¬ 
ing properties of signal transmission with time reversal. 

I. INTRODUCTION 

In a basic time reversal experiment, the intended receiver 
first broadcasts a short pilot pulse. The transmitter estimates 
the channel impulse response and sends the time reversed 
version of it back into the channel. These “time reversed” 
waves back propagate in the channel by retracing their paths 
and focus in space and time at the source, the intended 
receiver. 

In a channel with rich scattering, multipath or multiple 
scattering is exploited by TR [1] to focus broadband signals 
tightly in space and time. Moreover, the effective channel 
obtained through TR is hardened or statistically stable [2], 
Spatial focusing means that the spatial profile of the power 
peaks at the intended receiver and decays rapidly away from 
the receiver. Temporal focusing means that the channel im¬ 
pulse response at the receiver has a very short effective length. 
Channel hardening is the phenomenon of tightening of the 
distribution function of the effective channel impulse response. 
TR works well in systems where the delay spread of the 
channel times the bandwidth of the system is large [3], [2]. In 
this paper, we demonstrate TR with an Ultra-Wideband system 
in an indoor environment. 

The use of TR has three main benefits in communications. 
Temporal focusing significantly shortens the effective length of 
the channel. For example, the complexity of a MLSE equalizer 
is exponential in the length of the channel. TR, thus, reduces 
the complexity of the equalization task. 

A more important advantage of the TR technique is spatial 
focusing [2]. Spatial focusing results in very low co-channel 
interference in a multi-cell system. This results in a very 
efficient use of bandwidth in the overall system. 

Another advantage of using TR is the hardening of the 
effective channel, which means that TR results in a high 
diversity gain. The statistics of the time reversed channel 
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are different from the actual channel. Specifically the time 
reversed channel has a much smaller variance than the physical 
channel itself [2]. 

Fink and his coworkers [1], [3] (and references therein) 
have conducted an extensive number of experiments with TR 
in ultrasound and have shown large gains in spatial focusing 
in the presence of rich scattering. Also, TR methods have 
been used for underwater acoustic and ultrasound communi¬ 
cations [4], [5], [6], [3]. TR experiments using radio waves 
rather than acoustic waves have not yet been done. Here, 
we shall perform an initial demonstration of the effects of 
TR methods applied to radio wireless measurements taken 
from an indoor environment in the band 2-8GHz. Based on 
these measurements, we predict the performance of TR in 
such a setting. We demonstrate temporal and spatial focusing 
for a SISO link and evaluate channel hardening in terms of 
the empirical distribution of the received signal strength. Our 
prediction of the TR performance is exact if there is no channel 
measurement error at the transmitter and any imperfections in 
the transmit and the receive chain of the system are not too 
severe. 

In Section II we explain how data from static channel 
measurements can be used to predict the performance of a 
TR system, and we define suitable performance measures. In 
Section HI we describe how the measurements were made. In 
Section IV we discuss the results and finally conclude with 
Section V. 

II. TR Formulation 

Consider a transmitter-receiver pair. In TR the transmitter 
uses the time reversed complex conjugate of the channel 
impulse response as the transmit prefilter. Denote the channel 
impulse response by h{ ro, r), where r 0 is the receiver location 
and r is the delay variable. If the transmitter uses h*( ro, —t) 
as the prefilter, the effective channel to any location r is thus 
given by the time reversed field : 

s(r, r) = h*(r 0 , -r) * h(r, r) (1) 

where * denotes convolution with respect to delay. Note that 
a convolution with a time-reversed signal is equivalent to a 
correlation. 

In order to study the effects of TR we measure the channel 
impulse response between the transmitter and a receiver. We 
repeat the measurement by holding the transmitter fixed and 
changing the position of the receiver over a square grid. We 
take one comer of the grid as the reference receiver (the 
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intended receiver) and use its impulse response as the transmit 
prefilter /i*(r 0 , —r). The time reversed field at any point on 
the square grid can now be computed based on the channel 
measurements and by using Eqn. 1. This is equivalent to 
having the transmitter perform the prefiltering with the TR 
filter. 

The power of the signal as a function of both space and 
delay is computed and investigated. As a metric for spatial 
focusing we define the following quantity: 

k(t) = max |s(r, r)| 2 (2) 

T 

This quantity represents the power of the strongest tap at a 
receiver located at r. We expect «(r) to peak at ro and to 
decay rapidly with increasing distance from the origin. 

III. Measurement Setup 

Measurements were conducted by Intel Corporation at off- 
peak hours to ensure channel stationarity. The environment is 
an office space (40m x 60m) with many cubicles. Measure¬ 
ments span the bandwidth 2-8 GHz with 3.75 MHz frequency 
resolution. From the data we estimate that the coherence 
bandwidth of the channel is 20 MHz Antennas are vertically 
polarized. The virtual grid on which the receiver is moved has 
a distance of A 0 /4 where A 0 is the wave length of the mid 
frequency of the measurements (5 GHz). 

The receive antenna is moved to a different location with 
a precise robotic positioner. At each antenna position, the 
channel is measured using a vector network analyser. The 
measurements are corrected to compensate for the system 
components (including cable, gain stages, and antennas). The 
height of the transmit antenna is about 2.5m and that of the 
receive antenna is lm above the floor. The channel impulse 
response was computed by taking the inverse DFT of mea¬ 
surement data. 

IV. Experimental Results 

We evaluated «(r) for II different scenarios. Here we 
present two cases. Case I is an essentially line-of-sight scenario 
with a separation of 3m. Case II has a separation of 1 lm and 
is a typical non-line-of-sight situation. The results for the two 
cases discussed here have the worst focusing and yet they are 
very similar in their outcomes with all the other scenarios. 

In the 3-D figures, the square grid spans a region of 
7 Aq x 7Ao. A 3D plot of «(r) of the two cases is shown 
in Fig. 1. We see that spatial focusing works fine in both 
scenarios. In neither of them, however, the peak is isotropic; 
both peaks have one direction in which they Ml off faster, 
and another, in which the decay is slower. The structure of 
the peak carries some information about the geometry of the 
environment, i.e., about directions which show faster and those 
which show slower decorrelation in space. 

We also observe that the signal power level is at least 10 dB 
lower at a distance of 7A 0 than its value at the receiver. 
This demonstrates the spatial decorrelation very well. We can 
conclude that channel impulse responses show inherent quasi¬ 
orthogonality in space, which can be used for interference 
suppression. In Fig. 2 we demonstrate the compression of the 



Distance in wavelength Distance in wavelength 

(a) Case I 



„ Distance In wavelength 

Distance In wavelength 

(b) Case II 

Fig. 1. Spatial Focusing k ( r). One shot spatial field realizations for a) the 
line-of-sight b) the non-line-of-sight scenario. 


received pulse in the time domain. Displayed are the channel 
impulse response magnitude for the LOS scenario (a), the 
respective time-reversed, i.e., compressed impulse response 
magnitude (b), the impulse response for the NLOS scenario 
(c), and the compressed one for this case (d). The effective 
length of the channel is significantly shortened. Still, for 
the LOS-case, comparably strong side lobes remain visible. 
The impulse response of the NLOS channel is much better 
compressed in time. There is some intersymbol interference in 
the time compressed field and one would need to equalize the 
received signal. However, the time compressed channel has a 
significantly shorter effective length and is thus less expensive 
to equalize. 

In Fig. 3 we investigate the hardening of the channel. We 
consider the distribution of the strongest tap of the channel 
impulse response and of time-reversed field. The dashed curve 
represents the distribution of the zero-th tap of the channel 
impulse response, and the solid curve represents that of the 
compressed signal. The figure illustrates that the time reversed 
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Fig. 3. Distribution of s(r. 0) (solid) and that of /i(r,0) 


field is much less random (i.e,, has more diversity) than the 
channel response itself. 

V. Conclusions 

Using the indoor measurements discussed above, TR results 
in strong spatial focusing and time compression. At about 7 
wavelengths from the target the signal power reduces by at 
least 10 dB from its value at the target. The channel impulse 
response is well shortened at the target as well. Channel 
hardening or diversity gain is also observed. In a practical 
setting the results may degrade if the transmitter does not 
perfectly know the channel or if the transmitter circuitry has 
nonlinearities. Our calculations with the indoor radio wireless 
measurements obtained by Intel suggest that the benefits of 
TR, which have been observed in ultrasound and underwater 
sound TR experiments, could also be gained in wireless 
communications. 
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Fig. 2. Magnitude of channel impulse responses (a and c) and the time 
compressed impulse response (b and d) 





